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INTRODUCTION 


The  problem  of  atmospheric  transmission  of  infrared  and  visible  radi- 
ation Is  of  considerable  Interest  from  both  a purely  scientific  and  a practical 
point  of  view.  The  study  of  radiative  transfer  in  an  atmosphere  involves  a 
large  number  of  physical  and  chemical  effects  and  their  proper  theoretical 
treatment  requires  an  increasing  degree  of  mathematical  sophistication.  The 
vast  amount  of  physical  data  necessary  for  a detailed  description  is  typified 
by  the  AFGL  line  atlas  which  forms  the  core  of  the  AFGL  HITRAN  Research  Pro- 
gram*1^. Uses  of  this  data  Include  laboratory  studies  of  absorption,  computing 
data  for  comparison  or  design  of  field  observations  of  the  atmosphere  as  well 
as  such  applications  as  the  design  of  remote  sensing  or  laser  communication 
systems. 

In  the  design  and  simulation  of  such  systems,  the  atmospheric  trans- 
mission In  a given  wave  number  Interval  can  be  a crucial  consideration.  For 
systems  with  relatively  coarse  spectral  resolution,  the  detailed  spectral 
structure  of  the  absorption  Is  not  required.  Thus,  for  such  systems  a high 
resolution  technique  for  study  and  simulation  may  not  be  necessary  and  one 
may  profitably  use  lower  resolution  codes  such  as  the  AFGL  LOWTRAN  Program^, 
but  even  In  applied  work  a high  degree  of  spectral  resolution  is  increasingly 
in  demand. 

A number  of  line-by-line  calculatlonal  methods  have  been  reported 
which  provide  as  detailed  a spectral  resolution  as  computationally  possible. 
Among  these  may  be  found  the  LBL  Program  which  has  come  to  be  know  as  the 
HITRAN  Code^.  This  method  uses  the  AFGL  line  atlas  and  convolves  all  lines 
contributing  to  a given  wave  number  to  within  a prescribed  wave  number  range 
using  a Lorentz  line  shape  profile.  The  large  number  of  calculations  which 
the  LBL  Program  performs  requires  a large  amount  of  computer  time  which  makes 
parametric  studies  as  well  as  wlde-band  computations  prohibitively  expensive. 
Thus  this  code  tends  to  be  used  for  very  necessary  high  resolution  work  and 
also  as  a tool  to  improve  lower  resolution  models  such  as  LOWTRAN. 

With  the  growing  demand  for  high  resolution  work,  an  effort  was 
undertaken  at  AFGL  to  improve  the  efficiency  of  line-by-line  calculations. 

Tne  initial  phase  of  this  work  resulted  In  the  development  of  the  HIRACC 
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algorithm  for  the  convolution  of  Lorentz  line  shapes  within  48  half-widths 
of  a given  wave  number  for  a constant  atmospheric  temperature,  density  and 
composition.  This  extremely  efficient  algorithm  resulted  In  a decreased 
running  time  of  somewhat  more  than  an  order  of  magnitude  as  compared  to  the 
corresponding  parts  of  the  LBL  Code.  The  HIRACC  algorithm  has  been  fully 
documented  In  an  earlier  report^3-*  and  we  shall  assume  In  the  current  report 
that  the  reader  Is  familiar  with  Reference  [3]. 

At  this  point  It  would  be  helpful  to  the  reader  to  summarize  the  main 
features  of  the  HIRACC  algorithm  which  combine  to  give  effectively  an  order 
of  magnitude  decrease  In  computational  time.  The  first  step  In  the  develop- 
ment of  HIRACC  involved  the  decomposition  of  the  line  shape  profile  Into  three 
sub-prof des  with  different  spectral  widths.  The  algorithm  was  developed  for 
a pure  Lorentz  line  shape  and  the  line  shape  profile  was  cut  off  beyond  48 
halt-widths  from  the  line  center.  FASCODE  (Fast  Atmospheric  Signature  Code) , 
however,  was  designed  to  use  a Voigt  profile.  In  order  to  accomplish  this  it 
was  found  convenient  to  extend  the  cut-off  of  the  profiles  to  64  half-widths. 
This  enabled  the  definition  of  only  one  new  sub-profile,  a Gaussian  function 
representing  the  Doppler  contribution  out  to  4 half-widths.  Note  that  beyond 
this  number  of  half-widths  the  Doppler  broadening  may  be  neglected.  The  Voigt 
profile  is  then  approximated  by  an  appropriate  combination  of  four  functions. 
Figures  1.1  and  1.2  illustrate  the  resulting  decomposition. 

This  decomposition  has  the  advantage  that  the  rapid  variation  of  the 
profile  near  the  line  center  can  be  represented  properly  by  a specified  sampl- 
ing interval  and  the  more  slowly  varying  portion  as  one  approaches  the  wings 
of  the  line  can  be  properly  represented  with  a larger  increment  between  the 
points.  In  addition,  the  results  of  the  four  separate  convolutions  are  not 
put  together  until  the  contributions  from  all  the  lines  in  the  region  are 
complete.  An  interpolation  scheme  was  developed  to  compute  those  values  re- 
quired between  the  numerical  values.  Details  are  given  in  Section  2. 

Another  Important  result  of  Reference  [3]  was  the  determination  of 
a criterion  for  the  optimum  sampling  interval.  A sampling  interval  which 
gives  results  to  ...  0.1  percent  was  shown  to  be  given  by 

Av  - a/4  (1*1) 
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LORENTZ  DECOMPOSITION 


Figure  1.2:  Four  functions  used  to  reconstruct  the  Voigt  function. 


where  a Is  the  average  half-width  of  a line  over  the  region  of  Interest.  It 
was  also  shown  that  the  gain  G obtained  by  the  decomposition  method  over  com- 
puting the  entire  line  shape  at  a given  sampling  Interval  is  given  by  the 
relation 


G 


j nh  (sysx,,,) 

m*i 


(1.2) 


In  this  expression,  is  the  total  number  of  half-widths  over  which  the 
profile  is  required,  Nm  is  the  number  of  half-widths  spanned  by  the  mtfl  de- 
composed function  and  (6x./6x  ) is  the  ratio  of  the  sampling  interval  for  the 
first  function  to  that  required  for  the  nrn  function.  For  the  decomposition 
utilized  in  this  report,  this  yields  a theoretical  gain  of  5 1/3. 

We  also  point  out  another  feature,  namely  that  for  a given  spectral 
interval  (vj  - v2)  which  contains  N lines,  the  same  number  of  operations 
per  line  are  performed  for  the  sampling  Interval  chosen  according  to  the 
criterion.  This  Implies  that  the  convolution  of  a spectral  absorptance  for 
a given  spectral  Interval  will  always  take  the  same  amount  of  computational 
time  per  line  for  high  altitude  layers  where  the  lines  are  narrow  as  It  does 
for  lower  layers  where  the  lines  may  be  broadened  considerably.  This  feature 
is  borne  out  by  our  results. 

The  HIRACC  algorithm  has  been  used  to  develop  a multilayer  transmis- 
sion and  radiance  code  which  has  been  called  FASCODE  standing  for  Fast 
Atmospheric  Signature  Code.  In  the  course  of  producing  this  code  the  basic 
HIRACC  algorithm  has  been  left  intact,  but  several  peripheral  modifications 
have  been  implemented.  In  order  to  provide  the  capability  of  calculating  at 
higher  altitudes  where  the  Lorentz  profile  is  no  longer  appropriate,  a Voigt 
line  shape  version  has  been  developed  with  the  additional  feature  that  calcu- 
lations are  performed  out  to  64  half-widths.  One  of  the  great  advantages  of 
the  HIRACC  technique  is  t!  at  the  line  profile  may  be  changed  readily  with  only 
a small  change  in  running  time.  The  Voigt  line  shape  profile  is  described 
in  Section  2 of  this  report.  We  note  In  passing  that  purely  Lorentz  or  purely 
Doppler  versions  of  FASCODE  may  easily  be  implemented  for  the  user  whose  appli- 
cation is  in  these  domains.  We  emphasize,  however,  that  no  real  penalty  in 
running  time  will  be  paid  by  a user  who  exercises  the  Voigt  version  in  either 
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of  these  altitude  domains. 

As  mentioned  above,  the  HIRACC  algorithm  for  absorption  was  developed 
for  a uniform  atmospheric  path  (temperature,  pressure  and  absorber  concentra- 
tion). The  application  of  this  method  to  the  real  atmosphere  was  made  by 
approximating  the  real  atmosphere  by  a series  of  layers  with  constant  param- 
eters in  each  layer.  The  results  for  each  layer  are  then  merged  with  the 
absorption  coefficient  obtained  up  to  that  layer.  The  merging  algorithms  are 
described  in  detail  in  Section  3 of  this  report. 

Assuming  local  thermodynamic  equilibrium,  the  absorption  of  radiation 
by  a given  volume  of  gas  implies  the  re-emission  of  an  equal  amount  of  energy. 
Using  this  idea  one  may  "invert"  a transmission  calculation  to  obtain  the 
radiance  from  the  gas  along  the  given  path.  This  type  of  calculation  has  been 
provided  as  an  option  in  fASCODF.  to  enable  the  user  to  calculate  the  radiance 
due  to  the  atmsophere  itself.  For  the  user  who  has  a boundary  at  one  end  of 
his  path,  a provision  has  been  made  to  allow  the  user  to  add  his  own  surface 
radiance  model.  Currently  this  takes  the  form  of  a black  body  radiating  at 
a given  temperature  at  the  low  end  of  a path  starting  in  space  and  looking 
down.  This  could  be  easily  modified  if  a different  spectral  radiance  is  re- 
quired. This  part  of  FASCOOE  is  described  in  Section  4. 

In  Section  5,  we  present  some  avenues  for  future  efforts  and  summarize 
the  status  of  this  effort.  A series  of  appendices  gives  an  overview  of  the 
FASCODE  structure  as  well  as  complete  documentation  of  the  new  routines  in- 
volved. One  of  these  appendices  provides  a User's  Manual  for  some  sample 
inputs  and  outputs. 
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2.0 


LINE  SHAPE  - VOIGT  PROFILE 


To  calculate  spectral  absorption  contours  fo»  the  types  of  absorb- 
ing paths  encountered  In  the  atmosphere,  extensive  consideration  must  be 
given  to  the  spectral  line  shape.  For  atmospheric  paths  at  lower  altitudes, 
the  spectral  line  shape  within  ~5  cm"1  of  the  line  center  Is  given  satis- 
factorily by  the  pressure  broadened  (Lorentz)  shape.  Beyond  -5  cm"1  it  has 
been  shown  that  the  shape  1 r.  dependent  on  molecular  type  as  has  been  dis- 
cussed by  Burch**1,  Winters  et  a!*51  and  Clough  et  al*31.  At  high  altitudes 
(low  pressure)  the  line  shape  is  due  to  thermal  motion  resulting  in  the 
Doppler  line  profile.  The  a'.citude  above  which  the  Doppler  shape  is  valid 
is  dependent  on  the  wave  number  of  the  transition,  the  temperature  and  the 
molecular  species.  Between  the  domains  for  which  the  Lorentz  shape  and  the 
Doppler  shape  are  valid,  the  Voigt  line  shape  must  be  utilized.  Although 
the  Voigt  profile  was  originally  derived  for  the  case  where  broadening  was 
due  to  natural  broadening  and  thermal  motion  broadening,  it  is  applicable 
in  the  present  situation,  since  the  shape  function  is  the  same  for  natural 
as  for  pressure  broadening.  This  development  uses  an  approximate  method 
for  the  Voigt  profile  that  is  sufficiently  accurate  for  most  atmospheric 
problems  and  requires  substantially  reduced  computer  time  compared  to  pre- 
viously described  methods. 

The  absorption  coefficient  for  a pressure  broadened  line,  AL(v), 
as  a function  of  the  wave  number  value  of  the  radiation  field,  v,  the 
transition  wave  number,  vQ,  intensity  S and  half-width,  aL,  (half-width 
at  half  maximum:  HWHM)  is  given  by  the  Lorentz  function. 


A 


L 


qL 

+ (v  - v/ 


(2.1) 


The  pressure  broadened  half-width,  a,  , is  a function  of  the  absorbing  and 
broadening  molecule  types  and  the  vibratibnal -rotational  states  involved 
in  the  transition.  The  theory  that  has  proven  useful  for  calculating 
pressure  broadened  half-widths  is  due  to  Anderson^61,  as  implemented  by 
Tsao  and  Curnutte^71.  The  variation  of  the  half-width  with  temperature 
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and  pressure  Is  given  by  the  relationship 

\ <>•«  • \ <v  V ({)(;?) T (8-!) 

The  half-widths  on  the  AFGl  line  compilation  are  for  PQ  • 1 atm,  TQ  ■ ?98*K 
and  are  typically  of  the  order  of  0.08  cm’Vatm.  The  simplest  theoretical 
result  for  the  temperature  dependence  of  tls  half-width  gives  Xy  • 0.5. 

Recent  calculations  and  experiments  indicate  that  a higher  value  of  Xy  is 
more  appropriate,  Xy  ~ 0.75  (Benedict) The  absorption  coefficient  for 
a line  broadened  by  thermal  motion  is  given  by  the  Gaussian  function  in 
terms  of  the  Doppler  width,  cy  (HWHM)  as 


(2.J) 


(2.4) 


For  the  half-width,  ay  to  be  given  in  cm’1  at  temperature  T(*K),  c is  the 
velocity  of  light  in  cm/sec,  k is  the  Boltimann  factor  in  erg/deg,  NQ  is 
Avogadro’s  number,  and  M is  the  molecular  weight  of  the  molecule  type  in 
gr. 

Both  functions  have  been  defined  such  that  the  functional  value  at 
one  half-width  from  the  line  center  is  one-half  the  function  value  at  the 
line  center,  that  is 


A (v  ♦ «)  . 

~*7v^r  • i 

It  is  also  Important  that  the  integral  of  both  functions  yields  the  line 
strength,  that  is 


« 
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S * / A (v)  dv 


(2.6) 


For  the  purposes  of  this  development,  and  In  order  to  utilize  the 
algorithm  developed  by  Clough  et  al^,  It  will  be  useful  to  define  a di- 
mensionless argument  for  the  line  profiles, 


v - v. 


(2.7) 


so  that  the  Lorentz  profile  Is  given  by 


1 c 1 V - V 

AI  (*)  * i ~ — - — y With  z = - 

1 v aL  l7? 


(2.8) 


and  the  Doppler  profile  Is  given  by 


*0  <*> 


exp  - (tn2)  zc  with  z = 


v - v. 


The  Integral  and  half-width  properties  take  on  the  following  form 


(2.9) 


A (z)  adz 


(2.10) 


A (2  ■ ?)  . i 

T7i“«“5T  7 


(2.11) 


Several  excellent  programs  have  been  written  to  calculate  the  Voigt 
line  profile  which  may  be  regarded  as  the  convolution  of  the  Lorentz  function 
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with  the  Gaussian  function  (J.H.  Pierluissi,  P.C.  Vanderwood  and  R.B. 
Gomez^;  S.R.  Drayson^10^  and  B.K.  Armstrong ^ 11  *) . For  the  application 
of  performing  line-by-line  calculations  for  a multi-layered  atmosphere 
Involving  large  numbers  of  spectral  lines,  these  methods  require  a pro- 
hibitive amount  of  computational  time  and  provide  more  accuracy  than  is 
generally  required  for  such  problems.  Approximate  ns  for  the  computation 
of  the  Voigt  profile  have  been  suggested  by  Whiting*12^  and  Kielkopf*12^. 
The  method  described  in  the  latter  two  papers  Involves  approximating  the 
Voigt  function  by  a weighted  sum  of  the  Doppler  and  Lorentz  functions.  It 
is  an  extension  of  this  approach  that  we  have  utilized  in  this  development. 

The  star  ting  point  for  our  approximation  to  the  Voigt  function  Is 
the  definition  of  a Voigt  parameter  that  is  simply  calculable  and  is  well 
behaved  In  the  Lorentz  and  Doppler  limits.  For  a spectral  transition  with 
a Doppler  width,  aQ,  and  a Lorentz  width,  aL*  we  define  a Voigt  parameter, 
Ci  where 


C 


aL 

aL  + aD 


(2.12) 


In  the  Lorentz  and  Doppler  limits  respectively  we  have, 

C ■ 0 aL  <<:  a0  (Doppler  Limit)  (2.13) 


and 


C = 1 » aD  (Lorentz  Limit)  (2.14) 

The  definition  of  a yolgt  half-width,  ay,  is  now  required  to  proceed  with 
our  program.  This  definition  is  made  consistent  with  the  definition  of  the 
half-width  for  the  Dcppler  and  Lorentz  functions;  it  is  the  half-w  dt.i  at 
half  maximum  of  the  Voigt  function.  An  excellent  approximation  has  been 
given  by  Kielkopf  in  terms  of  the  Lorentz  and  Doppler  widths. 
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al 

o^  - (l  + e)  + 


0.0990  itn2. 


(2.15) 


In  order  to  utilize  this  expression  In  terms  of  c»  It  Is  necessary 
to  consider  the  determination  of  In  two  domains  of  c: 

o»v  * aVD^  * ao  0.0  < c < 0.5  (Doppler  Regime)  (2.16) 

and 


\ “ aVL 


U) 


a, 


0.5  < c <_  1.0  (Lorentz  Regime) 


(2.17) 


The  quantities,  aVD(c),  and  «yL(c).  may  be  obtained  from  Equation  (2.15) 
In  terms  of  c as 


end 


0.  < C <0.5 


11/2 

+ 1 

(2.18) 


«VL(c) 


0.5  < c < 1 . 0 


(2.19) 


The  largest  error  In  the  determination  of  the  Voigt  half-width  using  Equation 
(2.15)  Is  - 0.02%.  Values  of  avo  with  0.0  < c < 0.5  and  for  ayL(c)  with 
0.5  < c < K using  Equations  (2.18)  and  (2.19)  are  determined  for  equally 
spaced  values  of  c separated  by  0.005.  Plots  of  ayp  and  a^^  as  a function 
C appear  In  Figure  2.1  and  Figure  2.2  respectively.  This  error  is  not  as  great 
as  that  incurred  In  using  values  of  aVD  md  aVL  at  discrete  values  of  c for 
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values  of  ; Intermediate  to  those  tabulated.  For  most  physical  problems, 
the  errors  involved  in  the  present  method  of  determining  ay,  are  less  than 
the  errors  in  tne  values  of  and  aL  themselves.  If  it  is  desired  to  have 
more  accurate  values  of  ctyp  and  cty^,  the  Voigt  program  of  Armstrong  can  be 
utilized  in  an  Iterative  procedure.  In  general,  only  one  iteration  Is  re- 
quired to  give  a high  degree  of  accuracy. 

The  Voigt  function  may  now  be  approximated  as  a weighted  sum  of 
Doppler  and  Lorentz  functions  of  wldtn,  ay, 

AV(C.  V Z)  - (1  - Ck))  AD(av»  Z)  + C(c)  AL  (ay,  Z)  (2.20) 

with  Z ■ (v  - vQ)/ay 


This  definition  preserves  the  integral  property  d ~ired:  that  the  integral 
over  Z gives  back  the  line  strength,  S,  since  C(c)  is  independent  of  Z and 
the  two  functions,  AD  and  AL,  are  themselves  normalized  to  the  line  strength. 
Consider 

(X  OB 

f AVU,  v Z)  dZ  Ml  - C(c))  J A„(v  Z)  av  dZ 


and 


(ay,  Z)  ay  dZ, 


(2.21) 


* (1  - C(c))  S + C(C)S 


(2.22) 


so  that 


- S (2.23) 

The  remaining  problem  is  the  determination  of  the  weighting  constant,  C(c) * 

In  References  [12]  and  [13]  analytic  functions  have  been  given  for 
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the  constant  C(c)  as  a function  of  Voigt  parameter.  The  difficulty  with 
the  pro  edure  outlined  t.*ere  Is  the  V'.rgc  number  of  operations  required  to 
attain  C(c)  from  the  analytic  expressions  given.  The  most  straightforward 
way  to  proceed  Is  to  determine  the  numerical  function  C(t)  for  the  same  201 
values  of  t,  used  to  determine  the  Voigt  width.  Rather  than  using  the 
analytic  expressions  given  In  References  (12)  and  |U1  a least  squares 
procedure  has  been  used  to  determine  C(c)  by  minimising  the  weighted  sum 
of  the  deviations  squared,  a , as  obtained  from  Equation  (2.20), 


(2.24) 


All  the  quantities  have  been  defined  except  the  weighting  function  and 
the  grid  of  points  represented  by  2^.  Some  experience  Indicated  that  the 
weighting  scheme  which  presented  the  best  compromise  between  error  in  the 
central  portion  of  the  line  profile  and  the  line  wing  was  obtained  by  set- 
ting the  weight  to  the  inverse  of  the  Voigt  value,  that  is 


W1  " l/Ay(C»  tty.  Zj) 


(2.25) 


A more  obvious  choice  would  have  weighted  the  deviations  inversely  as  the 
square  of  the  Voigt  value  to  maintain  nearly  constant  percent  error  across 
the  line  profile.  It  was  deemed  more  Important  to  maintain  a smaller  per- 
cent error  for  the  larger  values  of  the  function  near  the  center  of  the 
line,  The  points  were  chosen  at  equally  spaced  Intervals  over  three  dif- 
ferent domains  for  reasons  that  will  be  discussed  later.  The  grid  was 
chosen  as  follows: 
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1-1-19  ?!  - 0,  Z2  = 0.25,  Z1S  ■ 4.0 

1 ■ 20,  30  ZgQ  * 5.0,  Zgj  ■ 6.0,  , Z^q  * 16.0 

1 " 31,  41  Z3j  ■ 20.,  Z32  “ 24.,  , Z ^ ■ 64.0 

The  results  obtained  for  ttu  constant,  C(c),  are  shown  In  Figure  2.3  as  a 
function  of  c-  The  values  are  determined  at  the  same  201  values  of  c In 
Figures  2.4  through  2.8.  The  continuous  curves  are  the  result  of  the  Arm- 
strong Voigt  program  and  the  x's  are  the  values  obtained  from  the  least 
squares  procedure.  For  c * 0 and  ; = 1 the  Voigt  function  Is  exactly  repro- 
duced as  a required  result  of  the  method.  The  results  are  given  on  a 
logarithmic  plot  in  order  to  give  perspective  to  the  value  of  the  function 
for  which  the  percent  errors  are  the  largest.  For  c = 0.05  to  c * 0.3,  the 
percent  error  In  the  wing  Is  of  the  order  of  22  percent,  but  these  errors 
occur  for  very  small  values  of  the  function.  Another  region  of  moderate 
percent  error  Is  at  Z ~ 3 half-widths  for  this  same  range  of  zeta.  The 
largest  percent  error  In  this  domain  Is  - 8%  for  zeta  - 0.2. 

As  Indicated  earlier,  these  results  are  sufficiently  accurate  for 
most  atmospheric  applications.  As  suggested  by  both  Whiting^12*  and  Klel- 
kopf^13^,  an  error  function  can  be  added  to  reduce  this  remaining  discrepancy 
to  negligible  proportions.  Some  discussion  of  this  will  appear  in  a subse- 
quent section. 

The  method  for  approximating  the  Voigt  profile  has  been  oeveloped 
a id  this  development  must  now  be  Incorporated  Into  the  algorithm  for  comput- 
ing spectral  absorption  coefficients  outlined  by  Clough  et  al^.  The  Voigt 
function  has  been  obtained  as  a weighted  sum  of  the  Doppler  function  and 
the  Lorentz  function  . For  purposes  of  convolving  the  approximate  Voigt  pro- 
file with  the  sf  ectral  line  data,  we  consider  Equations  (2.8)  and  (2.9). 

In  Reference  (3]  the  slow  convergence  of  A^(z)  has  bee.',  discussed  and  a 
method  was  oeveloped  to  reduce  computational  effort  for  performing  the  con- 
volition.  Fortunately,  the  Doppler  shaoe  has  very  rapid  convergence  to  zero 
as  a function  of  the  argument,  z The  domain  of  the  argument,  z,  for  which 
values  of  the  function  A0(i)  need  to  be  considered  are  limited  to  |z|  less 
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Figure  2.5:  Voigt  Line  Shape  Profile  for  t « 0.050 
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than  four  (|z|  £ 4),  noting  that  AQ  (z  = 4)  * 7.7  x 10'®.  The  normalized 
function  for  the  Doppler  line  profile  Is  given  by 


It  then  is  appropriate  to  decompose  the  Lorentz  function  A^(z)  into 
three  domains  of  z, 


0 < z £ 4 

0 < z < 16 
and 

0 < z < 64  (2.28) 


This  decomposition  Is  performed  in  a manner  similar  to  that  described  In 
’ Reference  [3].  A quartic  function  Is  defined  over  the  domain  0 < z < 4 
where 


Q,  = i («!  + b,  z2  ♦ c,  z4) 


(2.29) 


The  normalized  Lorentz  function,  L,  Is 


L-I  — y • 

71  1 + z2 


(2.30) 


The  coefficients  of  the  quartic  are  chosen  such  that  for  the  function  (L-  Qj) 
the  value  and  the  first  and  the  second  derivative  are  zero  at  the  boundary, 
Z^.  For  the  first  function,  7^  = Zj  = 4.  These  constraints  are  achieved  by 
the  following  relations 


aj-  (1  + 3 l[  + 3 zj)/(l  + Zj;) 


bj-  -0  + 3 Z^)/(l  + Z*)' 


(2.31) 

(2.32) 
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c « l/d  + Z*)3  (2.33) 


Similarly,  a second  quartlc  function,  Qp  Is  defined  such  that  for  function 
(L  - Q2),  the  value,  the  first  and  second  derivative  are  zero  at  the  bound 


ary  ■ Z2<  In  this  case  Zb 

the  function  utilized  Is  the 

obvious  in  tabular  form: 

* Z2  a 16.  In  the  region 
Lorentz  function  itself.  T 

from  16  £ |Z|  <6^ 
he  procedure  is  mor» 

DOMAIN 

Function 

0 £ |Z|  . 

< 4 

(Doppler) 

0 < |Z|  < 16 

0 < |Z|  < 64 

Xq(Z) 

G(z) 

(Lorentz) 

0 

0 

XL1<*> 

L(z)  - 

Qj(z) 

0 

0 

XL2(z> 

Qj(z)  - 

Q2(z) 

L(z)  - Q2(z) 

0 

XL3(z) 

Q2(z) 

L(z) 

; Note  that  the  Doppler  function  spans  the  same  domain  as  the  first  decomposed 

[ Lorentz  function,  that  the  functions  XL^  sum  to  L(z)  In  each  domain  and 

f that  the  functions  are  continuous  across  the  domain  boundaries. 

| The  four  functions  that  will  be  utilized  to  reconstruct  the  Voigt 

r function  are  shown  In  Figure  1.2.  The  functions  are  tabulated  at  201  values 

l of  the  argument  |z|  over  the  domain  valid  for  each  function. 

The  total  expression  for  approximating  the  absorption  due  to  the 
■ Voigt  profile  is  given  by 
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Avk.  z)  = A-j  (l  - C(C))  XD  (z) 
+ c(c)  [^Xu(z)  + XL2(z)  + XL3(z) 


(2.34) 


This  may  be  separated  into  three  functions  where 


Ha* 

II 

> 

3 

** 

k* 

> 

(l  - C(rj)  XD(z)  + C(c)  Xu  (z ) 

I 0 - 1*1 

< 4, 

iv2(s.  v z)  * 

c(0  xL2  (z) 

o i Ul 

< 16, 

iv3(C,  v z)  * -- 

c(c)  XL3(z) 

o 5 Ul 

*3- 

VC 

v| 

(2.35) 


The  sampling  interval  established  in  Reference  [3]  i dicates  that  Ayi  is 

aV 

sampled  at  intervals  of  -J-,  Ay2  at  intervals  ay,  and  Ay3  at  intervals  of 

4 *ay.  This  sampling  scheme  results  in  each  function  being  sampled  at  33 
values.  The  results  of  the  convolutions  over  the  three  domains  are  stored 
in  arrays  FF,  SF  and  VSF.  The  saving  in  computational  effort  results  from 
two  principal  reasons.  Only  99  values  are  required  to  describe  the  Voigt, 
function  over  + 64  half-widths.  The  composite  spectrum  from  the  three  arrays 
(FF,  SF,  and  VSF)  is  constructed  only  after  the  convolutions  have  been  com- 
pleted for  all  the  lines  in  a given  spectral  interval.  This  latter  p^nt 
is  one  that  was  not  sufficiently  stressed  in  Reference  [3].  The  implica- 
tion of  this  technique  is  that  not  only  is  the  grid  of  the  function  fine 
in  the  region  where  the  function  is  va-ying  rapidly  and  coarse  in  the  region 
where  the  function  is  varying  slowly  (actually  there  are  three  discrete 
sampling  intervals),  but  that  the  values  for  intermediate  points  are  not 
calculated  until  after  the  contributions  from  all  the  spectral  lines  have 
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3.0  TRANSMISSION  THROUGH  THE  REAL  LAYERED  ATMOSPHERE 

The  HIRACC  algorithm^  as  originally  reported  by  Clough  et  al  com- 
puted absorption  coefficients  for  a given  path  through  a gas.  The  thermo- 
dynamic properties  over  this  path  were  assumed  to  be  const  nt  which  enabled 
the  definition  of  a constant  sampling  Interval  over  the  spectral  region  of 
Interest.  The  path  Itself  was  defined  through  Input  In  the  form  of  column 
densities  of  the  appropriate  molecular  absorbers  together  with  the  temper- 
ature and  pressure.  We  note  that  column  densities  are  the  required  units 
for  the  AFGL  tape. 

At  first  glance,  It  might  appear  that  calculations  of  the  spectral 
absorptance  (or  optical  depth)  in  the  real  atmosphere  would  require  only  the 
definition  of  the  proper  values  of  these  column  densities.  However,  one 
must  take  Into  account  the  variation  of  thermodynamic  properties  along  the 
path  as  well  as  the  fact  that  certain  absorber  molecules,  such  as  H^O  and  03 
are  not  uniformly  mixed  at  all  altitudes.  For  example,  ozone  concentration 
peaks  In  the  stratosphere.  When  a calculation  of  the  optical  depth  at  a 
higher  altitude  is  required,  one  must  take  Into  account  the  fact  that  the  de- 
creased pressure  implies  a narrower  line  width.  Indeed,  as  one  goes  from  sea 
level  toward  space,  the  line  shape  profile  passes  from  almost  pure  Lorentz  to 
pure  Doppler.  If  one  uses  an  improper  sampling  Interval  to  calculate  the 
spectral  absorptance.  It  Is  clear  that  Important  absorption  features  can 
easily  be  neglected,  thus  yielding  Incorrect  results.  Therefore  a proper 
calculation  of  absorption  over  a path  traversing  a large  region  of  the  atmo- 
sphere will  require  an  appropriate  variation  of  the  sampling  interval  with 
altitude.  One  could,  of  course,  use  the  small  sampling  interval  required  at 
high  altitude  layers  for  all  layers,  but  this  would  require  unnecessary  cal- 
culations at  the  lower  altitudes. 

In  order  to  use  the  HIRACC  algorithm  with  as  little  modification  as 
possible,  we  have  taken  the  approach  discussed  by  McClatchey  et  al^l  and 
approximated  the  real  atmosphere  by  a series  of  layers,  each  defined  to  have 
constant  pressure  and  temperature  and  appropriate  values  of  the  column  den- 
sities of  the  absorbing  molecules.  Clearly  such  a decomposition  Is  not  unique 
and  requires  the  user  to  exercise  some  care  In  defining  the  atmosphere  as  will 
be  discussed  below.  While  this  arrangement  may  demand  a certain  degree  of 
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sophJstication  of  the  user,  it  Is  judged  th  . ,iis  solution  will  enable  a user 

to  be  more  flexible  In  changing  from  one  r-r<  .an  to  another  than  he  would  be 
if  the  program  were  tied  to  a given  atmospheric  model. 

The  basic  structure  of  FASCODE  is  then  a successive  application  of 
the  HIRACC  algorithm  for  successive  atmospheric  layers  with  an  appropriate 
merging  of  the  results.  For  convenience,  the  HIRACC  Program  reported  in 
Reference  [3]  has  been  left  almost  unchanged  except  to  modify  it  to  subrou- 
tine form  and  to  include  the  Voigt  line  shape  profile  described  in  the  pre- 
vious section. 

The  merger  of  the  abscrptance  for  two  different  layers  is  performed 
in  the  following  manner.  The  results  for  the  first  layer  are  computed  and 
written  to  disk  as  described  in  Reference  [3].  (We  remind  the  reader  that 
the  HIRACC  algorithm  processes  the  spectral  absorptance  in  "panels",  namely 
in  groups  of  data  defined  at  wave  numbers  separated  by  a sampling  interval. 
Typically  there  are  2400  quantities  in  a panel  but  a panel  may  be  shorter 
at  the  end  and  beginning  of  the  requested  wave  number  interval.)  The  reso- 
lution at  which  the  first  layer  is  to  be  calculated  is  determined  by  the 
sampling  interval  criterion  described  by  Clough  et  al^3^,  namely  one-quarter 
of  the  average  half-width  of  the  lines.  The  next  layer  to  be  computed  may 
have  pressure  and  temperature  sufficiently  different  such  that  a new  sampl- 
ing interval  is  required  by  the  sampling  criterion.  Let  the  sampling  inter- 
val for  the  first  layer  be  DV1  and  that  for  the  second,  DV2.  In  the  initial 
development  of  FASCODE  a decision  was  made  to  limit  the  program  to  specified 
ratios  DV1/DV2.  Note  that  since  the  pressure  in  the  atmosphere  is  a mono- 
tonically  decreasing  function  of  altitude,  the  average  half-width  of  a 
t*  ^al  Lorentz  line  profile  will  decrease  as  the  calculation  proceeds  from 
tc.  ower  to  higher  altitudes.  Thus  the  sampling  interval  determined  using 

the  /'iterion  of  Reference  [3]  will  also  decrease  as  the  calculation  proceeds 
to  the  higher  altitudes,  until  the  pressure  no  longer  remains  the  principal 
determ  ling  factor.  This  will  occur  at  altitudes  where  the  Doppler  width 
become  dominant  in  the  determination  of  the  half-widths.  Su'h  effects  are 
alreaoy  included  in  the  Voigt  line  shape  profile  model  discussed  above.  In 
addition,  the  Voigt  model  provides  a proper  treatment  of  the  transition  region 
at  those  altitudes  where  neither  the  Doppler  nor  the  Lorentz  profile  is 
appropriate. 
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In  general  then,  the  sampling  Interval  Is  larger  at  lower  altitudes 
and  smaller  at  the  higher  levels.  This  implies  that  a calculation  which 
proceeds  from  the  lower  to  the  higher  layers  Is  more  efficient  since  one 
would  not  have  to  perform  the  calculation  of  the  lower  layers  at  the  small 
sampling  intervals  required  by  the  narrower  line  shapes  at  the  higher  alti- 
tudes. Therefore,  in  FASCOOE,  the  HIRACC  algorithm  is  implemented  by  start- 
ing at  the  lowest  layer  (highest  pressure)  and  proceeding  to  higher  layers 
(lower  pressures). 

Since  the  decision  was  made  to  limit  the  number  of  ratios  of  the  sampl- 
ing intervals  between  two  adjacent  layers,  FASCODE  resets  the  DV  of  a new 
layer  to  the  nearest  allowed  ratio.  Note  that  the  user  must  be  careful  to 
prepare  the  atmospheric  input  with  a sufficient  number  of  layers  so  that  the 
criterion  for  the  ratios  is  met.  This  is  not  a difficult  condition  to  meet 
and  it  is  discussed  in  detail  in  the  User's  Manual  in  Appendix  B.  In  principle, 
the  restriction  to  a fixed  number  of  ratios  is  not  necessary.  This  was  dis- 
covered late  in  the  development  of  FASCOOE,  but  was  left  for  future  effort. 

We  shall  discuss  this  briefly  in  the  final  section. 

The  computation  of  the  second  layer  then  proceeds  using  the  HIRACC 
algorithm  with  the  reset  value  of  the  sampling  interval.  The  new  results 
are  also  written  to  disk,  by  Subroutine  PANEL,  using  a different  file  name. 
Having  obtained  the  spectral  absorptance  for  each  layer,  one  now  needs  a method 
to  merge  the  two  results  such  that  after  the  merger  one  has  the  absorptance 
for  a path  through  both  layers  defined  at  the  resolution  of  the  higher  layer. 
This  is  obtained  by  interpolating  the  "old"  or  coarser  resolution  tesults 
into  the  "new"  or  tiner  resolution  values  of  the  spectral  absorptance.  A 
simple  four-point  Lagrangian  interpolation  scheme  proved  adequate^1^ . This 
is  summarized  as  follows:  If  f(x)  is  a function  defined  numerically  over  a 
given  region  with  a constant  increment  h,  then  the  value  of  the  function  ac 
some  point  (xQ  + ph),  where  p is  a real  number  (-1  p < 2),  is  approximated 
by  the  relation 

f(xQ  + ph)  » A_jf_j  + AQf0  + Ajfj  + A2f?  (3.1) 

Here  f is  the  value  of  the  function  at  (x  + ph)  and  the  constants  are  given 
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by  the  relations 

A-i  * 

“P(P 

• 1)  (P 

- 2)/5 

(3.2) 

Ao  ■ 

<P2- 

1)  (P  - 

2)/2 

(3.3) 

A1  ■ 

-P(P 

♦ 1)  (P 

- 2 )/2 

(3.4) 

A2  ’ 

P(P2 

- 1 )/6 

(3.5) 

The  number  of  allowed  resolution  ratios  was  chosen  in  the  following 
manner.  First  it  was  required  that  the  old  ana  new  absorptances  should  be 
aligned  at  some  running  values  of  the  wave  number  at  least.  This  implied 
that  the  ratios  were  to  be  ratios  of  Inteqers.  For  further  simplification 
it  was  also  required  that  the  ratio  be  of  the  form  (N  + 1)/N.  In  addition, 
one  must  include  the  case  where  the  sampling  interval  does  not  change,  namely 
a one-to-one  ratio.  After  some  experimenting  with  larger  values  of  the  in- 
teger N,  it  was  found  adequate  to  limit  the  number  of  ratios  to  the  following; 
2/1,  3/2,  4/3,  5/4  and  1/1. 

The  interpolation  is  then  performed  by  identifying  the  quantity  xQ 
in  Equation  (3.1)  with  the  next  lowest  wave  number  of  the  old  array  below 
the  value  needed  for  addition  to  the  new  array.  The  interpolation  is  char- 
acterized by  an  index  (called  ITYPE  in  the  program)  which  is  the  number  of 
points  needing  interpolation  between  the  wave  numbers  which  co-align  in  the 
two  arrays.  This  will  be  clear  after  a glance  at  Table  3.1.  For  the  1/1 
case  (straight  add)  no  interpolation  is  required  and  the  index  is  zero.  For 
the  2/1  ratio,  one  point  must  be  Interpolated  and  ITYPE=1,  The  remaining 
schemes  are  quite  clear.  The  value  of  the  index  ITYPE  is  used  to  determine 
the  values  of  p in  Equations  (3.2)  - (3.5)  and  then  in  turn  to  compute 
arrays  of  the  interpolation  constants  for  interpolating  the  values  between 
the  co-aligned  wave  numbers.  The  remainder  of  the  merging  algorithm  consists 
in  bookkeeping  to  access  the  two  disk  files  containing  the  panel  data  for 
the  old  and  the  new  layers.  In  addition,  the  merged  results  are  also  written 
to  a third  disk  file  for  lerging  with  the  next  layer.  The  procedure  is 
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TABLE  3.1.  ALLOWED  RATIOS  AMD  SCHEMES 


then  continued  until  the  final  layer  is  processed. 

In  order  to  test  the  merging  algorithm,  the  following  procedure 
was  devised.  The  calculation  was  performed  as  described  above  except  that 
the  spectral  absorptance  calculated  for  each  layer  was  replaced  by  a known 
function  of  simple  form.  Consider  the  case  where  the  absorptances  from  two 
layers  are  merged.  If  vQ  Is  the  lower  boundary  of  the  frequency  range  of 
Interest  and  Av^and  Av2  are  the  sampling  Intervals  for  layers  1 and  2 re- 
spectively, (Avj  > Av2  by  assumption)  then  the  test  Is  made  by  replacing 
the  calculated  optical  depths  in  each  layer  (e.a  , Xj  and  x2)  by  the 
expressions 


Tj  ■ vQ  + (jj  - 1)  AVj/2  and  t2  = + (j2  - l)Av2/2  (3.6) 

where  jj  and  j2  are  running  indices  which  give  the  number  of  frequencies 
at  which  calculations  have  been  performed,  i.e.,  the  lower  bound  wave  num- 
ber is  j * 1 and  the  upper  bound  is  j = n^  where  n&  is  the  total  number  of 
frequency  values  in  the  region  calculated. 

The  array  is  to  be  merged  into  the  finer  spaced  array  x2.  Call 
the  merged  result  xJ2  and  let  j be  its  running  index.  From  Equation  (3.6), 
it  follows  readily  that 


|tj2  (j  + 1)  - x^2  (j)|  = Av2 


(3.7) 


provided  the  merging  process  is  correct.  Since  the  interpolation  involved 
cannot  be  expected  to  be  exact,  the  following  test  was  made  for  all  values 


I xi2  ^ 


t12  (j  - 1)| 


\^2  1 


10' 


(3.8) 


If  this  inequality  was  not  satisfied  the  program  was  directed  to  print  this 
fact  together  with  the  associated  parameters. 

This  testing  procedure  was  followed  and  the  inequality  (Equation 
(3.8))  was  found  to  be  satisfied  at  all  points  except  at  the  lower  and  upper 


38 


NV'^fVv. 


bounds  of  the  wave  number  interval.  This  was  due  to  the  problems  of  starting 
and  stopping  the  Lagrangian  interpolation  when  points  are  not  available.  A 
simplified  interpolation  scheme  was  chosen  at  these  boundaries  which  did  not 
have  sufficient  accuracy  to  satisfy  the  merging  criterion  in  Equation  (3.8). 
This  Is  not  a problem  however,  since  the  HIRACC  algorithm^  automatically 
expands  the  requested  wave  number  region  by  a small  amount  in  order  to  assure 
that  all  lines  which  contribute  to  a wave  number  interval  are  Included.  Thus 
we  conclude  that  the  merging  procedure  does  pass  our  criterion.  The  extension 
of  this  test  to  more  than  two  layers  is  straightforward.  In  Equation  (3.6), 
the  factor  1/2  is  to  be  replaced  by  1/N  where  N is  the  number  of  layers. 

The  next  test  of  the  algorithm  required  the  calculation  of  a case 

for  which  another  high-resolution  calculation  has  been  performed.  For  this 
purpose,  the  calculations  done  by  Kyle  at  NCAR  were  selected^ lf>^ . Kyle  per- 
formed a multilayer  atmospheric  transmission  calculation  in  the  wave  number 
region  (1-2600)  cm"*.  He  used  the  AFGL  tape,  a Voigt  line  shape  truncated 
at  5 cm"*  from  the  line  center  and  a model  atmosphere  based  on  the  AFGL  mid- 
latitude summer  profile^1**.  The  atmosphere  he  used  is  given  in  Table  3.2. 

We  note  briefly  that  Kyle's  results  were  degraded  in  resolution  by  convolu- 
tion with  a triangular  instrument  scanning  function  with  full  width  at  half 
maximum  of  20,  5,  and  0.1  cm-*.  Since  FASCODE  did  not  as  yet  have  this  fea- 
ture, we  compared  directly  with  his  highest  resolution  graphs.  It  is  also 

noted  that  Kyle's  atmosphere  had  to  be  modified  to  fit  the  ratio  criterion 

mentioned  above.  This  was  done  by  defining  additional  layers  where  there 
were  larger  altitude  gaps  in  Kyle's  atmosphere.  We  show  the  resulting  atmo- 
spheric model  in  Table  3.3.  We  limited  ourselves  to  the  wave  number  range 
(2000,  2200)  cnf*.  This  shall  be  referred  to  as  the  test  problem. 

The  results  of  the  calculation  for  our  test  problem  are  shown  in 
Figure  4.2.  The  spectra  obtained  by  FASCODE  are  in  entire  agreement  with 
those  obtained  by  Kyle.  We  delay  further  discussion  of  the  results  until 
the  next  section  in  order  to  presvent  transmittance  and  radiance  results  at 
the  same  time. 
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TABLE  3.2.  ATMOSPHERIC  LAYERS  USED  IN  THE  COMPUTATIONS  OF  KYLE  I16 3 
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TABLE  3.3.  ATMOSPHERIC  LAYERS  USED  IN  THE  TEST  PROBLEM 
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4.0  RADIANCE  FROM  AN  ATMOSPHERE  IN  LOCAL  THERMODYNAMIC  EQUILIBRIUM 

If  one  assumes  that  a given  infinitesimal  volume  of  a gas  is  in 
local  thermodynamic  equilibrium  (LTE)  at  some  temperature,  it  follows  that 
whatever  amount  of  radiant  energy  Is  absorbed  by  this  gas,  an  equal  amount 
of  energy  must  be  re-emitted  In  order  to  maintain  the  equilibrium  state. 

Furthermore,  the  spectrum  of  the  radiation  re-emitted  will  be  determined  by  ! 

the  black  body  or  Planck  function  using  Klrchhoff's  law.  This  LTE  model  may  ! « 

be  expected  to  be  more  nearly  valid  at  the  lower  altitudes  where  collision 
frequencies  are  quite  high.  At  the  higher  altitudes  (~  25  km),  one  should  j 

expect  non- equilibrium  effects  (NLTE)  to  be  Important.  Indeed,  Degges*17*  1 

has  been  developing  a comprehensive  NLTE  atmospheric  radiance  model  for  a 
number  of  years  and  other  workers  have  studied  the  high  altitude  problem  as 
well*-18-*.  The  FASCODE  Program  could  also  be  used  for  the  NLTE  case  if  vibra- 
tional or  rotational  temperatures  and/or  populations  of  states  at  higher 
altitudes  were  read  in.  Such  data  could  be  prepared  by  programs  such  as 
Degges'  work. 

Despite  its  shortcomings  at  higher  altitudes,  a LTE  atmospheric 
radiance  modil  can  be  very  useful  not  only  in  the  low  altitude  regime,  where 
it  rests  on  solid  footing,  but  also  in  the  higher  regions  where  one  may  use 
it  to  characterize  the  degree  of  deviation  of  the  NTLE  situation  from  the  j 

equilibrium  case.  A number  of  workers  have  prepared  LTE  radiance  models  among 
which  we  want  to  mention  the  recent  extension  of  the  AFGL  LOWTRAN  model  to 
include  a radiance  calculation^. 

In  this  section  we  present  the  algorithms  implemented  in  FASCODE  to 
enable  the  calculation  of  atmospheric  radiance  along  a given  optical  path 
assuming  LTE  along  that  path.  To  provide  the  most  efficient  calculation  of 
the  radiance,  it  was  decided  to  calculate  its  value  layer-by-layer  along  with 
the  spectral  absorptance  calculation  described  in  the  previous  section.  It 
should  be  noted  that  the  radiance  of  a path  proceeding  from  lower  toward 
higher  levels  is  not  the  same  as  that  for  the  same  geometrical  path  proceed- 
ing from  higher  to  the  lower  layers.  Figure  4.1  presents  a sketch  of  two 
such  paths,  spacw-to-ground  (Figure  4.1a)  and  ground-to-space  (Figure  4.1b). 

A four-layer  atmosphere  has  been  shown  for  simplicity,  each  layer  being 
labeled  by  the  letter  A,  B,  C,  or  D.  The  boundaries  between  each  layer  have 


Figure  4.1:  Sketch  of  the  difference  between  a path  looking  up  and  one 
looking  down 
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been  numbered  from  1 to  5 where  5 Is  the  ground  at  temperature  at  0^.  The 
temperature  of  each  layer  Is  called  0^  ^ where  i Is  one  of  the  labels  for 
the  boundaries.  The  transmittance  contributions  between  any  two  boundaries 
are  defined  as  T^  where  1 and  j are  boundary  Indices. 

Consider  now  the  equilibrium  radiance  from  Layer  A to  space  In  Fig- 
ure 4.1a.  Using  Klrchhoff's  law  this  may  be  written  as 

RA-$  * (1  ’ T12}  P(012’  v)  (4.1) 

*■  i 

where  P Is  the  Planck  function  and  v Is  the  frequency  In  cm  . Proceeding 
to  the  next  layer  the  contribution  of  this  layer  to  the  radiance  observed  In 
space  Is 


B-3 


12 


,(i  - t23) 


P ( ®23  * ^ 


and  the  rema  nlng  two  layers  can  be  written  as 
RC-$  * T12  T23  ^ " T34^  P^034* 

" T13^  " T45^  P^034’  v) 


(4.2) 


(4.3) 


and 


T„  ^9a(l  “ ^45)  r^®45» 


'D-S  '12  '23  '34 


(4.4) 


* T14{1  ‘ T45>  P(045*  v) 


Finally  the  contribution  of  the  ground  Is 


Rg-s  " 1 12  T23  T34  T45  P^0g’  v)  * T15  P'9g* 


(4.5) 


The  total  radiance  Is  then  the  sum  of  all  these  .erms.  Extending  this  to 
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the  case  of  n layers  with  n + 1 boundaries  we  can  write 


R(down)  „ 
total 


n 

A 


1,1 


(1  - T 


1.  1+1 


me, 

I f 


+ T 


1’  n+1  rvvn+l 


P(e«*i.v) 


(4.6) 


where  now  en+1  Is  a temperature  characteristic  of  the  last  boundary.  Note 
that  in  general,  this  is  not  necessarily  the  ground.  (It  shall  be  noted 

J 

rij  a ^ Tk  k+r * A1so*  1^  Is  not  required  that  the  spectral  distribution 

at  the  boundary  be  given  by  the  Planck  function.  An  arbitrary  spectrum  or 
one  characteristic  of  a particular  type  of  radiating  boundary  may  be  substi- 
tuted for  the  (n+1)' st  term. 

Turning  now  to  the  other  case,  looking  up,  we  can  write  the  follow- 
ing expression  for  the  contribution  to  the  radiance  from  each  of  the  layers 


R5-4  " ^ " 745^  P^045’  v' 

(4.7) 

R4-3  = T45^  " T34^  P^034’ 

(4.8) 

R3-2  = T34  T4S  ^ " T23^  P^023’ 

(4.9) 

“ T35  ^ “ T23^  p(023’  ^ 


R2-l  * T23  T34  T45  ^ ' TlP  P^012’ 


"1*25  ^ ~ 


(4.10) 
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Again,  the  total  radiance  Is  the  sum  of  the  Individual  terms  for  each  layer 
which  may  be  written  as 

M "I1  , x . x 

"(total)  = Jo  T"-1'  » (!  ' T"-'-1-  »-<)  P(en-1-1.  "-t  V) 

(4.11) 

+ Tl.  n+1  p<el-  v) 

where  we  have  Implicitly  defined  ^ = 1 for  convenience.  To  present  a radi- 
ating boundary  at  the  end  of  an  upward  looking  path,  such  as  a cloud  or  the 
Zodiacal  light,  one  may  add  the  same  type  of  boundary  radiating  term  as  dis- 
cussed for  the  downlooking  case.  An  upward  looking  boundary  has  not  been 
Included  In  FASCODE  as  yet  but  a user  can  readily  add  one  if  it  should  be 

required.  A glance  at  Equations  (4.6)  and  (4.11)  is  sufficient  to  see  the 

difference  between  each  case. 

The  expressions  in  (4.6)  and  (4.1J)  are  convenient  mathematical  repre- 
sentations of  the  algorithm  but  for  computational  purposes  it  is  more  useful 
to  represent  the  algorithms  such  that  they  manipulate  the  new  increment  to  the 
transmittance  for  a given  layer  i,  and  the  radiance  and  total  transmittance 
accumulated  up  to  that  layer.  As  was  mentioned  in  Section  3,  the  calculation 
proceeds  from  the  lowest  layer  to  the  highest  in  order  to  minimize  the  time 
for  merging.  If  AT.j  is  the  incremental  transmittance,  E(1)  the  radiance  and 
P(i)  the  Planck  function,  the  radiance  and  transmittance  after  the  i**1  layer 
has  been  traversed  are 

(up)  (up) 

E(1)  = E(i-l)  + (1  - ATi ) P(i)  • T(i-l)  (4.12) 

and 

(down)  (down) 

E(1)  = AT,  E(i-l)  + (1  - ATi)  P(i)  (4.13) 

for  the  radiance  looking  up  and  looking  down,  respectively.  Th°  transmittance 
is  clearly  given  by  the  relation  (for  both  cases) 
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T(1)  - AT1  • T(1-l) 


(4.14) 


The  boundary  radiance  is  frequently  not  a function  with  very  fine 
spectral  resolution.  Thus,  for  the  lookingdown  case,  the  most  efficient  pro- 
cedure Is  to  compute  the  boundary  radiance  at  the  coarse  resolution  of  the 
lowest  layer  and  Include  It  with  the  radiance  of  the  first  layer  (modified  by 
the  transmittance,  of  course).  This  will  then  be  properly  carried  along  through 
the  remainder  of  the  calculation.  For  the  case  of  a radiating  boundary  at  the 
upper  end  of  a path  looking  toward  space,  it  is  simplest  to  add  the  contribu- 
tion from  the  boundary  after  all  calculations  have  been  performed  modifying 
It  by  the  total  transmittance  computed  for  the  path.  As  mentioned  above, 

FASCODE  at  this  juncture  does  not  Include  coding  for  this  second  case,  but 
does  Include  a radiating  boundary  at  the  end  of  a down-looking  path. 

In  order  to  speed  up  the  calculation  of  the  Planck  function  the  fol- 
lowing procedure  was  followed.  For  each  panel,  the  black  body  function  is 
computed  and  the  results  stored  in  an  array  with  separation  of  one  wave  number. 
The  value  needed  at  a given  wave  number  is  obtained  by  Interpolation  which  is 
performed  only  when  the  value  of  the  Planck  function  can  be  expected  to  have 
changed  sufficiently  to  warrant  an  updated  value.  If  the  Planck  function  is 
written 


P(v)  » Av3  (esv  - 1) 


(4.15) 


where  v is  in  wave  numbers  and  s = C2/0»  0 being  the  temperature,  and  C2 
the  second  radiation  constant  (1.4388).  Taking  the  derivative  with  respect 
to  v we  may  write 


dP 

T 


dv 

V 


SV\-1 


3 - (sv)  (1  - e ) 


(4. If) 


If  one  takes  | AP/P | ~ 10"*,  the  increment  at  which  one  should  Interpolate  to 
obtain  a new  value  of  the  Planck  function  is  easily  computed.  This  alleviates 
the  need  to  recalculate  the  Planck  function  unecessarily. 
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The  Implementation  of  this  algorithm  for  LTE  radiance  In  FASCODE  Is 
contained  In  two  subroutines,  EMUP  and  EMDOWN,  for  paths  "looking"  up  toward 
pace  and  "looking"  down  toward  the  ground  respectively.  An  additional  sub- 
routine was  prepared  to  compute  the  first  layer  (the  lowest),  EMINIT.  In 
this  routine  the  possibility  of  radiation  from  a boundary  was  Included  for 
the  downward- looking  case. 

The  merging  of  the  results  layer-by-layer  proceeds  In  a similar 
fashion  to  that  for  the  absorptance  described  In  Section  3.  One  difference 
is  that  now  one  needs  to  merge  two  quantities,  namely,  the  radiance  computed 
and  the  transmittance  as  the  calculation  proceeds  from  layer  to  layer.  These 
quantities  are  written  to  disk  In  two  records,  first  the  radiance,  followed 
by  the  transmittance.  The  code  accesses  the  HIRACi.  algorithm  to  obtain  the 
spectral  absorptance  panel  by  panel  for  each  layer.  This  Information  Is  con- 
verted to  transmittance  by  exponentiation  in  either  EMUP  or  EMDOWN  depending 
on  the  case  of  interest.  The  LTE  radiance  algorithm  is  then  exercised  as  the 
merging  Is  taking  place. 

The  radiance  routines  were  tested  in  a manner  similar  to  the  merging 
test  described  In  the  last  section.  Consider  the  case  of  a series  of  n layers 
with  a path  from  space  to  ground  with  a boundary  at  temperature  0.  If  one 
modifies  the  input  data  such  that  each  layer  has  the  same  temperature  as  the 
boundary,  then  one  can  show  that  the  resulting  radiance  will  be  given  by  the 
Planck  function  at  temperature  0,  a result  which  is  not  at  all  surprising. 

The  same  result  holds  for  the  upward  looking  case  If  a boundary  Is  added  at 
the  upper  end  and  the  temperatures  are  all  set  equal.  This  procedure  was 
followed  with  the  results  as  expected  to  within  accuracy  requirements  Inherent 
In  tr,e  method. 

Thf  idlance  package  has  been  tested  for  tne  problem  described  in 
previous  section.  The  radiance  results  are  shown  in  Figures  4.2  together 
wltf;  the  transmittance  for  a path  looking  from  space  to  ground  which  Is 
modeled  as  a black  body  radiator  with  temperature,  T = 273  K.  Note  the  smaller 
range  of  the  abscissa  for  Figures  4.2(e)  and  4.2(f).  The  absorption  features 
seen  In  Figures  4.2  are  clearly  seen  in  the  radiance  profiles  as  well.  The 
expanded  frequency  scale  In  Figures  4.2(e)  and  4.2(f)  show  this  most  clearly. 
Self-reversal  '"n  ea,  be  seen  in  the  radiance  for  strong  absorption  lines. 
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Figure  4.2  'e):  Transmittance  and  radiance  for  the  test  problem  with  expanded  scale  (2060-2065  cm 


nee  and  radiance 


This  is  particularly  dramatic  in  Figure  4.2(f)  near  the  band  edge  of  the  COg 
(11101-00001)  transition. 

In  Table  4.1  we  give  the  time  results  of  the  run  made  to  produce 
Figures  4.2.  All  times  are  for  the  CDC  6600  computer  system  at  AFGL.  For 
each  of  the  layers  of  the  test  problem,  we  give  the  altitude  boundaries  and 
the  temperature  and  pressure  characterizing  each  layer.  The  next  column  gives 
the  value  of  ; used  In  computing  the  Voigt  shapes  for  that  layer.  The  average 
linewldth  and  the  sampling  Interval  are  giv^  in  the  columns  labeled  ALPHA  and 
DV  respectively.  The  ratios  of  the  sampling  intervals  are  given  In  the  follow- 
ing column.  The  final  four  columns  show  the  timing  for  the  convolution  of  the 
spectra  (C0NV),  thewrlting  of  the  panels  to  disk  (PANEL),  the  calculation  of 
the  radiance  (EM)  and  the  calculation  of  the  transmittance  merging  (ABS).  The 
sums  of  the  times  for  each  column  are  also  presented.  Note  that  the  total 
number  of  lines  processed  was  6681.  The  convolution  took  approximately  0.7 
msec  per  line  per  layer.  This  latter  statistic  is  very  meaningful,  since  it 
gives  the  reader  some  idea  of  the  speed  of  the  program.  The  total  time  for 
a transmittance  calculation  only,  may  be  obtained  by  summing  the  totals  of 
the  columns  marked  C0NV,  PANEL  and  ABS  (i.e.,  102.2  seconds).  For  the  radiance 
calculation  along  the  same  path,  one  obtains  the  total  time  by  summing  the 
totals  of  the  columns  marked  C0NV,  PANEL  and  EM  (i.e.,  175.3  seconds). 


TIMING  RESULTS  FOR  THE  TEST  PROBLEM 


5.0 


CONCLUSIONS  AND  RECO  FNDED  EXTENSIONS  OF  THE  PROGRAM 


For  certain  purposes  it  is  helpful  to  use  a pure  Lorentz  or  a pure 
Doppler  line  profile  instead  of  the  Voigt  profile,  which  takes  slightly  more 
computational  time  than  the  Lorentz  line  shape  or  the  Doppler  line  shape 
alone.  For  example,  a user  might  have  a problem  involving  only  transmission 
at  the  very  high  layers  of  the  atmosphere,  40  km)  in  which  region  the  line 
shape  Is  purely  Doppler.  On  the  other  hand,  one  might  be  interested  in 
studying  a laboratory  experiment  at  relatively  high  pressures  for  which  the 
Lorentz  profile  is  adequate.  The  addition  of  the  pure  Doppler  and  the  pure 
Lorentz  cases  is  straightforward  and  has  been  accomplished.  The  reader  is 
referred  to  the  program  listing  in  Appendix  3.  The  revisions  which  were  re- 
quired for  the  pure  Lorentz  case  include: 

a.  Subroutine  SHAPED  is  not  needed. 

b.  The  least  squares  fits  for  the  linear  com- 
bination of  the  two  profiles  are  not  neces- 
sary (AVRAT,  aVD(c),  ayL(c),  C(c)). 

c.  All  program  references  to  the  quantity  c 
are  deleted. 

d.  The  final  result  is,  of  course,  similar  to 
the  original  HIRACC  coding^.  Here  it  is 
called  HIRACL. 

For  the  pure  Doppler  case,  the  following  differences  arise: 

a.  Subroutine  SHAPEL  is  not  used. 

b.  Delete  all  c references. 

c.  Only  the  FF  array  is  used,  VSF  and  SF  are 
not  needed. 

d.  The  resulting  routine  is  called  HIRACD. 

Some  care  had  to  be  taken  to  assure  that  the  proper  indexing  is  made. 

The  contributions  of  the  various  continuum  features  have  not  been  in- 
cluded as  yet-1’2’1*’20^.  These  features  arise  from  a number  of  physical 
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processes  such  as  absorption  by  atmospheric  aerosols,  and  a variety  of  molec- 
ular continuum  processes  for  molecules  such  as  Ng,  COg,  H^O.  In  addition  the 
contributions  of  the  lines  beyond  64  half-widths  must  also  be  Included  as  part 
of  the  continuum  contribution.  In  general,  the  continuum  absorption  Is  a slow 
function  of  wave  number  and  may  be  directly  Incorporated  into  the  VSF  array. 
This  operation  will  not  significantly  affect  the  running  time  of  the  program. 

The  line-by-line  results  have  finer  spectral  detail  than  Is  required 
for  comparison  with  some  experiments.  Some  Instrumentation  (especially  that 
used  for  engineering  systems)  does  not  have  resolution  such  that  the  final 
spectral  detail  computed  in  this  version  of  FASCODE  can  be  resolved.  A method 
is  required  for  convolving  a given  instrument  scanning  function  with  the 
FASCODE  output  in  order  to  degrade  the  detailed  results  for  comparison  with 
lower  resolution  data.  Work  on  this  aspect  is  currently  underway  at  AFGL. 

We  note  that  the  convolution  techniques  used  in  the  line-by-line  spectral 
synthesis  (HIRACC  algorithm)  may  also  be  applied  to  a scanning  function 
convolution. 

When  the  sampling  interval  is  constant  from  one  layer  to  the  next, 
and  one  is  not  calculating  the  radiance,  it  seems  clearly  possible  to  devise 
a new  method  which  would  decrease  the  running  time  considerably.  This  method 
would  compute  at  the  same  time  the  spectrum  for  all  of  those  adjacent  layers 
which  have  a constant  sampling  interval.  For  example,  above  a certain  alti- 
tude where  the  Doppler  line  shape  becomes  dominant,  essentially  all  of  the 
sampling  intervals  can  be  taken  to  be  constant  and  the  merging  calculations 
for  all  of  these  layers  can  be  done  panel  by  panel. 

Finally,  for  applications  in  systems  studies  it  is  recommended  that 
spherical  geometry  be  added  which  could  account  for  the  fact  that  the  earth's 
atmosphere  is  not  plane  pa.allel,  but  rather  spherical.  With  this  addition 
to  the  code,  the  calculation  of  limb  radiance  would  be  easily  performed  and 
practical  applications  could  be  taken  directly  from  FASCODE. 

To  illustrate  this,  consider  a remote  sensing  satellite  viewing  the 
earth  on  some  sight  path.  If  the  sight  path  ends  on  the  earth's  surface  or 
some  other  surface  above  the  earth  such  as  a cloud  layer,  no  great  extension 
to  the  program  need  be  made.  One  merely  needs  to  program  a method  for  com- 
puting column  densities  of  the  absorbing  molecules  for  each  layer.  However, 
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for  the  case  of  a sight  path  that  Is  glancing,  or  tangent  to  some  altitude 
above  the  earth,  the  nroblem  is  more  complicated. 

Figure  5.1  is  a sketch  of  this  case.  Here  rmax  is  the  earth-centered 
radius  to  the  maximum  altitude  for  which  calculations  are  to  be  made,  and  rT 
is  the  corresponding  radius  to  the  tangent  height  at  Point  T.  Let  A and  3 be 
the  points  at  which  the  llne-of-slght  Intersects  the  circle  of  radius  rmax. 

First,  we  consider  the  case  where  only  the  transmittance  is  to  be  calculated. 

It  is  clear  that  the  largest  sampling  interval  which  will  be  determined  from 
the  criteria  programmed  In  FASCODE,  will  be  found  in  the  tangent  height  layer 
and  the  smallest  will  be  for  the  layers  ending  at  the  two  points  A and  B. 

Thus  the  optimum  procedure  will  be  to  start  at  the  tangent  height  where  one 
has  the  smallest  number  of  points  needed  to  characterize  the  convolved  spectral 
transmittance.  Next,  we  note  that  the  two  paths  ta  and  TB  are  identical  compo- 
sitions of  column  densities.  Thus  the  transmittance  from  path  TB  Is  identical 
to  that  on  path  TA  and  the  total  transmittance  may  be  obtained  simply  by  doubl- 
ing the  contributions  from  each  layer  and  computing  only  one  of  the  two  paths 
TA  or  TB.  Note  that  It  Is  the  column  densities  which  are  to  be  doubled  and  not 
the  transmittance.  Also  we  have  tacitly  assumed  that  the  atmospheric  composi- 
tion profiles  do  not  change  appreciably  along  the  path  ATB.  Since  the  angle 
ACB  can  be  as  large  as  - 20°,  this  may  not  be  the  case.  Thus,  for  example, 
the  sight  path  may  enter  the  atmosphere  at  - 45°N  latitude  and  exit  at  65*^ 
latitude  and  the  profiles  can  be  quite  different,  especially  for  water.  This 
Is  a complication  which  may  or  not  bo  Important  in  a given  case.  No  essential 
difficulties  should  occur,  however,  should  this  situation  need  to  be  investigated. 

For  the  radiance  calculation,  It  is  again  best  to  start  at  the  tangent 
part  of  the  sight  path.  This  can  be  done  simultaneously  with  the  TB  part  cor- 
responding to  a path  looking  up  and  the  TA  part  to  a down-looking  path.  When 
the  calculations  for  the  two  final  layers  are  finished,  the  total  radiance, 

Rjqj,  Is  computed  simply  by  the  relation 

RT0T  “ RTA  + TTA  RTB  (5,1) 

where  and  are  the  results  of  the  radiance  calculations  on  paths  TA 
and  TB  respectively,  and  T^  is  the  transmittance  from  path  TA.  If  Point  B 


61 


■ . — -n™;--  u ipy^n.ii.i^Hn  W)fw.  m,  HIM  i J.'iuu,  ■■>J1M  JUlH'i-gHUm 


Is  not  at  the  sanu  altitude  as  Point  A (for  example  Point  B corresponds  to 
a target  or  a cloud),  no  essential  problems  occur.  One  merely  stops  the 
calculation  along  TA  at  the  altitude  of  Point  B,  performs  the  composition 
of  Equation  (5.1)  and  continues  the  calculation  to  Point  A as  before,  but 
using  the  value  RTQT  for  the  radiance  at  the  stopping  point.  This  method 
has  not  been  Implemented  as  yet  In  FASCODE  but  the  Implementation  should  be 
straightforward. 
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In  this  apperdlx  we  present  a description  of  the  code.  For  the 
reader's  convenience  much  of  the  documentation  of  the  HIRACC  algorithm^ 
has  been  reproduced  here  with  an  occasional  modification.  In  addition,  the 
new  subroutines  developed  In  this  effort  are  described. 

The  program  has  been  written  to  use  line  parameter  input  data  con- 
sistent with  that  contained  In  the  AF6L  line  parameter  tape^.  The  line 
data  has  been  reformatted  onto  a binary  file  which  contains  the  line  data 
pertinent  to  the  molecules  and  wave  number  range  of  interest.  This  step 
has  been  taken  to  keep  read  time  consistent  with  the  time  required  to  per- 
form the  calculatlonal  part  of  the  program.  The  control  parameters  are 
read  from  the  input  file  and  written  to  the  output  file;  TAPE3  is  the  binary 
file  containing  the  line  parameter  data;  and  TAPE12  is  always  the  binary 
output  file.  The  final  output  and  Intermediate  output  files  contain  a 
header  record  which  includes  the  identification  information,  SECANT,  temper- 
ature, pressure,  molecular  identification,  and  molecular  column  densities 
of  the  homogeneous  layers.  The  first  record  for  each  output  panel  is  a 
header  record  for  the  pinel  which  contains  the  wave  number  values  of  the 
first  and  last  absorption  coefficient  values  of  the  panel,  the  wave  number 
increment  between  output  points  and  the  number  of  output  points.  The  second 
record  of  the  panel  contains  either  the  array  of  absorption  coefficient 
values  resulting  from  the  convolution  when  only  absorptance  is  calculated  or 
the  radiance  array  followed  immediately  by  the  transmittance  array  for  radi- 
ance calculations.  The  current  version  of  the  program  outputs  a maximum  of 
2400  values  per  panel;  in  general,  the  first  and  last  panel  are  shorter. 

The  output  file  MFILE,  and  an  additional  file,  LFILE,  are  used  to  store  the 
intermediate  data  prior  to  merging  the  results  for  the  layers.  The  merged 
results  are  always  on  MFILE.  The  file  KFILE  contains  the  absorption  co- 
efficients for  each  layer.  Table  A1  outlines  the  t.se  of  these  file*  and  file 
labels. 

FASCODE  has  been  designed  to  be  used  with  the  CDC  segmentation  fea- 
ture which  allows  the  program  to  be  run  using  a minimum  of  central  memory 
by  loading  dynamically  only  those  subprograms  which  are  in  use  at  a given 
stage  of  the  code's  execution.  Using  this  feature,  FASCODE  requires  only 
- 56  Kq  words  of  central  memory.  If  a user  does  not  have  a segmentation 
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capability  on  the  machine  being  used,  allowance  will  have  to  be  made  for 
a larger  amount  of  central  memory.  In  this  case,  storage  allocation  should 
be  redistributed  to  minimize  total  storage  requirement,  estimated  at  77  Kg. 
That  Is,  long  extension  arrays  in  EMUP,  EMDWN,  and  ABS  should  be  changed  to 
unlabelled  common-. 

A simple  overall  structure  of  the  code  was  obtained.  This  is  shown 
schematically  in  Figure  A.l.  HIRACV  and  Its  associated  subroutines  perform 
the  spectrum  synthesis  of  the  absorptance.  The  merging  subroutine  is  called 
ABSMRG  for  the  absorption  coefficient  case.  The  LTE  radiance  computations 
are  performed  by  subroutines  EMINIT,  EMUP,  and  EMDOWN.  Subroutine  TPLOT  is 
an  expanded  version  of  the  plotting  program  described  in  Reference  [3],  and 
it  has  been  changed  into  subroutine  form.  Note  that  FASCODE  can  be  used  to 
prepare  plots  directly.  Some  possible  future  extensions  of  the  program  have 
been  sketched  as  dotted  lines  in  Figure  A.l.  The  reader  should  note  that 
the  modular  construction  of  the  code  allows  straightforward  extension  and 
revisions. 

The  program  consists  of  the  main  Program  FASCODE:  Subroutines  ABSMRG, 
TPLOT,  EMINIT,  EMUP,  EMDOWN,  HIRACV,  SHAPED,  SHAPEL,  MOLEC,  RDFILE,  CONVFNV, 
PANEL,  HIRACL,  CONVFNL,  HIRACD.  CONVFNC,  and  PANELD;  and  the  Function  QVRFAC . 
The  overall  strategy  of  the  HIRACV  Subroutine  is  indicated  in  Figure  A. 2.  All 
the  subroutines  are  called  from  the  main  program  or  from  their  subdriver  (See 
Figure  A.l)  and  the  flow  of  the  program  is  easily  traced.  Subroutines  HIRACL 
and  HIRACD  perform  the  same  functions  for  the  Lorentz  and  Doppler  line  shape 
profiles  respectively. 

The  main  Program  FASCODE  reads  the  input  data  and  calls  each  of  the 
subdrivers  shown  in  Figure  A.l,  according  to  the  particular  run  desired  as 
determined  from  the  input.  After  initializing  constants  and  reading  the 
basic  parameters  for  the  run,  the  program  enters  a loop  for  the  calculation 
of  the  requested  results  for  each  layer.  Inside  the  loop  the  atmospheric 
properties,  average  temperature,  average  pressure  and  absorber  molecule 
column  densities  are  read  from  the  input  file  for  the  given  layer.  The 
proper  sampling  interval  is  then  computed  and  an  identification  header  for 
the  layer  is  buffered  out  to  KFILE.  Note  that  extensive  use  of  BUFFER  IN 
and  BUFFER  OUT  is  made  to  increase  efficiency.  If  a user  does  not  have  this 


70 


*mmmt **  wmjwuu  .♦ 


f ENTRY 

_n: 


Initialize  Program  Parameters 


T 


CALL  SHAPEL 
CALL  SHAPED 


CALL  MOLEC  (1) 
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Figure  A. 2:  Flow  Diagram  for  HIRACV  Subroutine 
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Test  to  Output  Panel 


CALL  PANEL 


Test  For  Calculation  Complete 


Test  To  Continue  Convolution 


SUBROUTINE  PANEL 

- Interpolates  VSF  Array 
Into  SF  Array;  SF  Array 
Into  FF  Array 

- Outputs  Panel  to  KKFILE 

- Resets  Arrays 

- Sets  Indicator  for 
Calculation  Complete 


Figure  A. 2:  Flow  Diagram  for  HIRACV  Subroutine  (Continued) 
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capability,  these  statements  can  bt  replaced  by  binary  read  ana  write  state- 
ments. Subroutine  NIRACV  Is  then  accessed  to  compute  the  spectral  absorptance. 
Once  this  has  been  completed,  the  appropriate  merging  routine  is  accessed, 
unless  It  Is  the  first  layer  for  which  no  merger  Is  necessary.  For  the  radi- 
ance calculation,  Subroutine  EMINIT  1$  called  to  compute  the  radiance  due  to 
the  first  layer.  For  all  other  layers,  the  radiance  Is  obtained  hy  calls  to 
EMDOWN  or  EMUP  for  the  cases  of  looking  do^n  and  looking  up  respectively. 

When  all  layers  have  been  processed,  the  plotting  Subroutine  TPLOT  is  called 
If  the  plot  flag  has  been  set  by  the  user.  Subroutines  SHAPEL  and  SHAPED  set 
up  the  convolution  functions  used  to  define  the  Voigt  function  from  0 to  64 
half-widths. 

Subroutine  KOLEC  In  conjunction  with  Function  QVRFAC,  makes  the 
molecular  Identifications  associated  with  the  line  parameter  file,  and  de- 
termines the  correction  factors  for  the  line  Intensities  (SCOP)  and  the  half- 
widths (ALFCOR).  The  quantity,  SCOR,  Is  the  correction  factor  due  to  the 
temperature  dependence  of  the  vibrational  and  rotational  partition  sums.  The 
vibrational  partition  sum  Is  calculation  for  a given  molecular  type  as 

N 

VT)  " 1-1  „-S7kTjdl  <A‘li 

where  ^ Is  a fundamental  frequency  and  d^  Is  the  degeneracy  of  the  vibration, 
The  temperature  dependence  of  the  rotational  partition  sum  Is  given  by 
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where  F ■ 1 for  linear  molecules  and  1,5  for  nonlinear  molecules.  The  refer- 
ence temperature,  T„,  Is  taken  as  296  K,  consistent  with  the  AFGL  Line  11st- 

o 1211 

Ing.  For  further  discussion  of  these  topics,  see  Herxberg'  pp.  503  ff. 

The  partition  sum  calculations  are  performed  by  QVRFAC  and  the  necessary 
molecular  parameters  are  contained  In  data  statements  In  Subroutine  MOLEC. 

The  quantity,  ALFCOR,  Is  the  correction  factor  due  to  the  pressure  and  temper- 
ature dependence  of  the  collision  broadened  half-width.  The  temperature 
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dependence  of  the  half-width  has  been  taken  as  (T0/T)0,5  although  calculations 
based  on  the  Anderson-Tsao-Curnutte  theory  are  reported  to  give  somewhat  dif- 
ferent temperature  dependencies  (Varanasi) *22^. 

Subroutine  RDFILE  reads  the  blocked  binary  line  parameters  over  the 
wave  number  range  for  which  line  data  Is  required.  The  line  parameters  in- 
clude the  wave  number  value  of  the  transition  (GHU,  cm-*),  the  intensity  of 
the  transition  at  296#K  (S,  cm-*/ (mol/cm2)),  the  collision  broadened  half- 
width at  half  maximum  for  296°K  and  1 atm  pressure  (Al.FAO,  cm-*),  the  lower 
state  energy  of  the  transition  (EPP,  cm-*),  and  the  molecule  identification 
number  (MOL).  If  the  line  parameter  data  is  insufficient  to  complete  the 
specified  calculation,  the  message  "end  of  file  on  disk"  is  printed  on  the 
output  file.  If  no  further  line  data  is  required,  IDATA  is  set  to  1,  and 
control  is  returned  to  the  main  program. 

At  this  stage  of  the  main  program,  an  effective  optical  depth  is  cal- 
culated for  each  line  which  is  dependent  on  the  column  density  of  the  layer, 
the  secant  of  the  angle  through  the  layer,  the  temperature  of  the  layer,  and 
the  half-width  of  the  line  oi(ALFI). 

The  effective  depth,  u' » 


where  w is  the  absorber  column  density,  E"  is  the  lower  state  energy,  Q^R  = 
QyQR>  and  the  othe10  quantities  have  been  previously  defined.  In  terms  of 
the  program  coding  the  effective  depth  appears  as: 


EFDPTH  = SEC*W*S*SCOR*P.ECALF* 

EXP(EPP/XKTFAC)  + (l.-EXP(-GNU/XKT))/(l. -EXP(-8NU/XKT0)) 

where  XKT  and  XKTO  are  the  wave  number  values  of  T and  TQ  in  cm-*.  As 
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h previously  discussed,  the  proper  sampling  Interval,  DV,  should  be  0.?5  times 

|!j  the  average  line  half-width.  If  the  half-width,  ALFi,  in  less  than  the 

[|  sampling  interval,  the  half  width  Is  set  to  the  sampling  Interval  and  a 

“ series  of  minus  signs  Is  written  to  the  output  file.  If  the  half-width 

b exceeds  a maximum  value  (ALFMAX)  where  ALFMAX»60UND/S4  and  BOUND  Is  the 

S;  • 

s'  maximum  value  In  wave  numbers  over  which  a line  can  be  calculated,  the 

». 

(I  half-width  is  reset  to  ALFMAX,  and  a series  of  + signs  is  written  to  the 

output  file.  Included  in  the  records  indicating  the  resetting  of  the  half- 
r width  is  the  wave  number  value  of  the  transition  (GNU),  the  intensity  (S) 

| ! and  the  half-width  (ALFAO)  values  of  the  transition  from  the  line  parameter 

(' , cape,  the  calculated  value  of  the  half-width  (ALFI),  the  value  to  which  the 

| half-width  has  been  reset  (DV  or  ALFMAX),  and  molecular  identification  num- 

ber (M).  If  the  number  of  ha'if-wldth  changes  (NCHNG)  exceeds  100,  the  com- 
< putation  is  terminated. 

| Subroutine  CONVFNV  Is  a tightly  written  subprogram  in  which  consider- 

i able  effort  has  bee;;  taken  to  minimise  operations  In  the  DO  30  loop.  This 

[ subroutine  perfosns  tho  triple  convolution  of  XF,  XS,  ana  XVS  with  a line 

datum  putting  the  results  In  the  proper  elements  of  FF,  SF,  and  VSF  respect- 
: Ively.  A simplified  flow  diagram  appears  in  Figure  A. 3.  Control  indicator 

J.PANEL  Is  set  to  IDATA  if  the  DO  loop  over  the  lines  (40)  is  satisfied  indi- 
cating whether  e panel  is  complete  or  more  lines  are  required.  If  the  line 
DO  loop  (40)  is  not  completed,  IPANEL  1$  set  to  1 indicating  that  a complete 
panel  has  been  calculated.  Control  is  returned  to  the  main  program. 

If  IPANEL  has  been  set  to  1,  Subroutine  PANEL  is  called.  Subroutine 

' t 

PANEL  performs  a four-point  Lagrange  interpolation  of  the  VSF  array  into  the 
SF  array  and  the  SF  array  into  the  FF  array,  thus  combining  the  results  of 
three  independent  convolutions  into  a final  result.  A general  flow  diagram 
of  PANEL  is  given  In  Figure  A. 4.  Care  Is  taken  to  store  array  values  required 
for  the  Interpolation  of  subsequent  panels.  VFT  Is  the  wave  number  value  of 
the  first  element  of  the  FF  array,  which  is  common  to  the  first  element  of 
the  SF  and  VSF  arrays.  A binary  header  record  is  written  to  the  binary  file 
(KFILE)  for  each  panel  which  includes  the  wave  number  value  (VIP)  of  the  first 
element  of  the  panel  (FF(NLO)),  the  wave  number  value  (V2P)  of  the  last  ele- 
ment of  the  panel  (FF(NHI)),  the  wave  number  increment  (DV),  and  the  number 
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Figure  A. 3:  Flow  Diagram  for  SUBROUTINE  CONVFNV. 
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Figure  A. 4:  Flow  Diagram  for  SUBROUTINE  PANEL. 
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of  absorption  coefficient  values  outputted  (NLIM).  The  second  binary  record 
contains  the  NLIM  values  of  the  absorption  coefficient  from  the  FF  array.  The 
arrays  are  appropriately  shifted  and  reset  in  preparation  of  the  computation 
of  subsequent  panels.  Control  Is  again  returned  to  the  main  Subroutine  HIRACV. 

At  the  conclusion  of  the  outputting  of  the  last  panel  for  each  layer, 
a record  is  written  to  the  output  file  Indicating  the  current  value  of  the  time, 
the  time  spent  In  RDFILE,  in  LONVFN,  and  in  PANEL  (the  units  are  seconds). 

Also  Included  In  this  record  are  the  first  and  last  wave  number  values  of  the 
panel.  A second  record  is  written  to  the  output  file  indicating  the  average 
value  of  the  half-width,  the  number  of  lines  read  since  the  last  panel  was 
completed,  and  the  total  number  of  lines  read  since  the  initiation  of  the  con- 
volution calculation.  Control  is  returned  to  Statement  10  if  the  calculation 
is  complete,  or  to  Statement  40  to  continue  the  convolution  in  process. 

ABSMRG  — Absorption  Merging 

Subroutine  ABSMRG  performs  the  merging  of  the  absorption  layer-by- 
layer  along  a path.  The  absorption  at  the  present  layer  and  the  accumulated 
absorption  from  previous  layers  are  read  from  disk  files.  The  constant  ISMAIL, 
computed  in  FASCODE,  is  used  to  identify  whether  the  present  or  the  accumu- 
lated layer  has  the  smaller  wave  number  increment.  BUFFER  IN  and  BUFFER  OUT 
are  used  rather  than  READ  and  WRITE  for  speed  of  running.  The  array  FF 
contains  the  absorption  from  the  smaller  increment  and  DUMF,  the  absorption 
from  the  larger.  The  index  ITYPE,  computed  in  FASCODE,  gives  the  value  of 
the  numerator  of  the  ratio  of  the  two  increments.  The  coefficients  for  the 
Lagrange  four-point  interpolation  are  computed  and  an  interpolation  is  made 
for  tho  absc-pticn  in  the  layer  with  the  larger  increment.  The  merging  of 
the  two  layers  is  then  carried  out  such  that  the  merged  absorptance  has  an 
increment  equal  to  the  smaller  sampling  interval.  The  output  file  contains 
a header  record  containing  identification  information,  the  values  of  the 
secant,  pressure  and  temperature  concentrations  of  the  absorbing  molecules, 
wave  number  increment,  the  value  of  the  first,  and  last  wave  numDer  in  the 
panel  and  a layer  count.  Each  panel's  output  consists  of  two  records.  The 
first  record  contains  the  first  and  last  wave  number  in  the  panel,  the  wave 
number  increment  and  the  number  of  points  in  the  panel.  The  second  record 
contains  the  accumulated  absorption.  For  speed  of  computation,  if  the  wave 
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number  increment  is  the  same  for  both  layers,  the  merging  is  done  directly 
by  merely  adding  the  corresponding  values,  A simplified  flow  diagram  is 
shown  in  Figure  A. 5. 

Subroutine  EMINIT  computes  the  LTE  radiance  for  the  first  altitude 
layer  only.  The  Planck  function  is  computed  as  described  in  the  text.  There 
is  also  the  option  of  the  radiance  from  a ooundary,  for  the  case  of  looking 
down.  A simplified  flow  cha^t  Is  given  in  Figure  A. 6.  The  radiance  for  the 
case  looking  from  the  space  to  ground  Is  computed  by  the  formula 

NEW EM  « (1  - TR)*BB  + TR*0LDEM 
The  radiance  from  ground  to  space  uses  the  formula 

NEWEM  * (1  - TR)*BB 

where  TR  is  the  incremental  transmission  of  the  layer,  BB  the  Planck  func- 
tion ano  OLDEM  the  radiance,  from  the  boundary  when  it  is  requested.  Other- 
wise, it  Is  zero.  The  result^  of  the  calculation  are  written  on  a file  con- 
sisting of  four  records: 

1.  A header  record  containing  alpha-numeric 
information  supplied  from  FASCODE  in  the 
fi^st  seven  words,  the  values  of  the  secant, 
c.  ,e  pressure,  average  temperature,  con- 
centrations of  the  absorbing  molecules,  DV, 
first  and  last  wave  numbers  and  a layer  count. 

2.  VI,  V2,  DV  and  the  number  of  points  in  the 
panel . 

3.  The  radiance. 

4.  The  transmission 


Arrange  File  Names;  Read  Headers 
From  New  and  Old  Files;  Write  to  Merge  File 


Straight 


Read  Data  on  Both  Files.  Panel  by 
Panel;  Add  and  Write  to  Merge  File 
Read  When  All  Panels  Completed 


Set  Up  Interpolation  Coefficients 
Read  Absorption  Data  From  New  File. 

Read  Absorption  Data  From  Previous  Merge  File. 
Perform  Interpolation;  Write  to  Merge  File. 
Continue  While  Testing  For  Need  to  Read 
In  New  Panel  From  Either  File. 

When  a New  Panel  Must  be  Read,  Save  Data  At 
End  of  Current  Panel  For  Proper  Interpolation; 
When  All  Panels  Done,  Return. 


Figure  A.  5;  Flovf  Diagram  for  SUBROUTINE  ABSMR6. 
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Read  Header  Record 


Read  Panel  Data 
Set  Up  Planck  Function  Array 


Radiating"' 

.Boundary/ 


Set  Up  Boundary 
Radiance 


Set  Up  Criterion  For  Interpolating  Planck  Array 
Read  First  Layer  Absorption  Data 
Compute  Transmission  by  Exponentiating 
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Compute  Radiance 
Looking  Up 


Compute  Radiance 
Looking  Down 
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Figure  A. 6:  Flow  Diagram  for  EMINIT 
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EMUP  — Radiance  and  Transmittance  Looking  Up 

Subroutine  EMUP  performs  the  merging  of  the  radiance  and  transmis- 
sion layer-by-layer  along  a path  looking  up  from  ground  to  space.  As  in 
Subroutine  ABSMR6,  the  1/1  ratio  of  the  DV's  (wave  number  increment)  Is 
done  separately  for  speed  of  computation.  For  the  remaining  ratios,  ITYPE 
is  computed  in  FASCODE  and  equals  the  value  of  the  numerator  of  the  ratio 
of  the  two  DV’s.  The  coefficients  for  the  Lagrange  four-point  interpolation 
are  computed  and  a file  is  read  which  contains  the  infoimatlon  for  the  accumu- 
lated quantities.  This  file  consists  of  a header  record  which  includes  Identi- 
fication information,  the  value  of  the  secant,  pressure  and  temperature, 
concentrations  of  the  absorbing  molecules,  wave  number  increment,  the  value 
of  the  first  and  last  wave  number  !n  the  layer,  and  a layer  count.  This  is 
followed  by  a series  of  three  records  per  panel.  The  first  record  contains 
the  first  and  last  wave  number  in  the  panel,  the  wave  number  Increment  and 
the  number  of  points  in  the  panel.  The  second  record  contains  the  accumulated 
radiance  which  is  stored  in  array  OLDEM.  The  third,  the  accumulated  transmis- 
sion stored  In  array  OLDTR. 

Similarly,  the  file  for  the  new  layer  contains  the  header  record 
and  first  and  last  wave  number  record.  The  next  record  contains  the  spec- 
tral absorptance  from  which  the  transmission  is  computed  and  is  stored  in 
array  TR. 

The  black  body  function  is  computed  as  described  in  the  text  and 
an  Interpolation  is  made  for  the  radiance  and  transmission  with  the  larger 
wave  number  increment.  The  total  radiance  along  the  path  length  is  computed 
using  the  formula 

NEWEM=OLDEM  + (1.0  - TRi )*BB.*0LDTR 

where  OLDEM  is  the  interpolated  old  radiance,  TR..  the  new  transmission  incre- 
ment, BB^ , the  black  body  and  OLDTR  the  interpolated  transmission.  The  trans- 
mission is  computed  as 

NEWTR=OLDTR*TR. 
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The  output  file  is  written  such  that  It  can  be  used  as  the  input  when  treat- 
ing the  next  layer.  The  flow  chart  for  Subroutine  EMUP  is  in  Figure  A.  7. 

EMDOWN  — Radiance  and  Transmittance  Looking  Down 

Subroutine  FMDOWN  performs  the  merging  of  the  radiance  and  transmis- 
sion layer-by-layer  along  a path  looking  down  from  space  toward  ground.  How- 
ever, to  achieve  efficiency  it  is  arranged  such  that  the  calculation  always 
starts  at  the  lowest  altitude  and  continues  upward.  The  logic  for  EMDOWN  is 
exactly  the  same  as  Subroutine  EMUP  and  the  formula  for  computing  transmission 
is  the  same.  The  LTE  radiance,  however,  is  computed  as  follows: 

NEWEM-(1  - TRi)*BB1  + TR^OLDEM 

where  OLDEM  is  the  interpolated  old  radiance,  TR^  the  transmission  increment 
of  the  i'th  layer  and  BB^  the  black  body  function.  In  looking  from  space  to 
ground  it  is  not  strictly  necessary  to  compute  the  total  transmission  in  order 
to  compute  the  radiance  but  this  Mature  has  been  included  in  this  subroutine 
to  make  it  compatible  with  Subroutine  EMUP.  Figure  A. 7 also  gives  the  tiow 
diagram  for  this  routine. 
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Set  Up  Lagrangian  Interpolation  Coefficients.  1 

Read  Radiance  and  Transmission  From  Previous  Merge  File.  } 


Read  Absorption  Data  for  New  Layer. 

Set  Up  Planck  Array  for  Interpolation. 

Exponentiate  to  Get  New  Transmission  Increment. 
Interpolate  Old  Radiance  and  Add  New  Contribution 
Using  Up(Down)  Radiance  Formula  . 

Interpolate  Old  Transmission  and  Multiply  by  New 
Transmission  Increment. 

Write  Results  to  Merge  File  For  Use  With  Next  Layer. 


NO  X/V 
■— < Panels 

Done> 


Figure  A. 7:  Flow  Diagram  EMUP 
EMDOWN 


RETURN 
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Appendix  B presents  samples  of  input  and  output  from  a run  of 
FASCODE  on  the  AFGL  CDC  6600  computer.  A user  may  find  these  pages  help- 
ful in  assessing  whether  the  program  has  been  adapted  properly  to  a new 
machine  environment. 

As  pointed  out  already,  extensive  use  has  been  made  of  the  BUFFER 
IN  and  BUFFER  OUT  capability  available  in  CDC  FORTRAN  Extended.  Since 
other  computer  systems  (such  as  IBM  360/370)  may  not  have  this  feature 
available,  we  have  collected  a list  of  all  the  BUFFER  statements  together 
with  their  identification  numbers.  This  Is  given  In  Table  B.l  and  should 
be  of  assistance  to  a user  whose  machine  requires  the  changing  of  the 
BUFFER  statements  to  binary  READ  and  WRITE  statements. 

Table  B.2  gives  a sample  job  setup  for  a run  used  in  developing 
Test  Case  Output.  The  first  record  contains  CDC  job  control  cards  (SCOPE) 
This  is  followed  by  a series  of  Segmentation  Loader  input  cards  which  de- 
termine the  way  in  which  FASCODE  is  to  be  loaded  into  the  system.  The 
third  and  final  record  contains  the  input  data  for  the  execution  of  the 

Test  Case  Problem  by  FASCODE.  Table  B.3  contains  a list  and  description 

1 

of  the  FASCODE  input  data.  A load  map  for  this  run  is  found  in  Table  B.4. 
Finally,  the  output  pri.ited  by  FASCODE  in  executing  this  job  is  presented 
in  Table  B.5.  To  ascertain  that  a faithful  version  is  in  hand,  it  is 
recommended  that  a user  execute  this  test  case  and  compare  the  output  with 
Table  B.5. 

With  regard  to  the  preparation  of  model  atmospheres  for  input  to 
FASCODE,  the  principal  requirement  Is  to  assume  that  ratios  of  the  sampl- 
ing intervals  between  a layer  and  the  one  just  above  it  is  less  than  or 
equal  to  2/1.  It  is  suggested  that  3/2  would  be  a good  target  ratio.  Thl 
can  be  checked  by  computing  the  average  half-width  of  the  layers  using 
Equations  (2.2)  or  (2.4).  If  a user  would  feel  more  comfortable  in  his 
mind  that  a good  model  has  been  derived,  one  or  two  extra  layers  could  be 
added  at  relatively  low  cost.  However,  it  should  be  clear  that  each  new 
layer  increases  the  running  time  correspondingly. 
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TABLE  B.l.  LOCATION  OF  BUFFER  STATEMENTS 


BUFFER 

IN ( L INF  XL , 1)  (FILHO(l) ,VL INHI) 

U C k 6 3 L 

BUFFER 

IN  (UNFIL.l)  (FILHr(t),  VLINMI) 

000790 

BUFFER 

OUT (KFIL  E , 1)  (710(1) ,Nl»YER> 

002540 

10 

SUFFER 

IN  (LINFIL,1>  <VHIN,NREC) 

01 3660 

BUFFER 

IN  (tlNFIU,!)  < C.NU(Ii>  ,GNU(I2)  ) 

0037 10 

3C 

BUFFER 

IN  (LINFIL,!)  ( 6NU(I1),M0KLIMIN>  > 

00J75t 

BUFFER 

OUT  (KFILi.l)  (VIP.NLIM) 

005250 

BUFFER 

OUT  (KFILE,!)  (FF (NLO) »FF (NHI ) ) 

315270 

8UFFER 

IN(KFILE,1)  ( X I f! ( 1 ),  XI 0(2)  ) 

0090  40 

260 

9UFFER 

OUTMFIUE  ,1)  (XIOU)  i NL4YE  9) 

0C9Q60 

420 

BUFFER 

IN (L FILE, 1) (XIO(l) ,NL*YER) 

CC9000 

BUFFER 

IN(IOLn,l) <OVlP,NLINO> 

0C  9190 

BUFFER 

XN( XOLO , 1 ) (OllMF (1)  , OUMF  (NL IMO - i 

C49210 

BUFFER 

IN  i.  KNEW,  1 ) (VIP,  NL  1*0 

0C  92  30 

SUFFER 

IN ( KNEW , 1 ) (FF (1) , PF(NL IN)  ) 

009250 

buffer 

OUT (MFILE, 1) ( V IP , NL IN) 

3C9290 

BUFFER 

OUT(MFTLE,l>  (FF (1 ) ,FF  (NL I M) ) 

0C 9310 

BUFFER 

TN(IOLB,l) (OV1 P(NLIHO) 

019550 

jin 

SUFFER 

INCOLB,  t)  (OUMF(l)  ,OUMF(NUMO>  ) 

009580 

BUFFER 

IN  ( KNEW  , 1 ) (V1P.NLIM) 

01  9620 

33  0 

BUFFER 

IN  (KNEW,  1 ) (FF(  1) , FR(NLIS)  ) 

UC965G 

BUFFER 

IN( IOLB , 1)  (O V 1P,N LINO) 

009920 

BUFFER 

IN ( TOLO , a ) (BUMF (4 ), OUMF (NL TNO > ) 

009970 

me 

BUFFER 

OUT(NFILE,l)  (V1P,NLIN) 

01C170 

370 

BUFFER 

OUT (NFlLE , 1) IFF(t) ,FF(NLIM>) 

010190 

mo 

buffer 

XU (KNEW, 1 ) (V1P,NLIM) 

C1025C 

750 

BUFFER 

IN(KNEN,1) (FF (1) ,FF (NL IN ) ) 

I)1C  20w 

buffer 

IN(KFILE, 1) (XI 0(1) , NL4YER) 

010720 

BUFFER 

OUT (NFILE ,1) (XIO(I) ,NL«YE») 

010740 

BUFFER 

IN (KF ILE , 1 > (OV1P.NLTHO) 

OlUdO 

BUFFER 

IN<  KFILE , 1)  ( FF (1)  , FF ( NL I M0) ) 

* 1 1«  5S 

buffer 

OUT (MFILE , l) (OVIP.NLIMO) 

011360 

BUFFER 

OUT  (MFIL  E,  1)  (ENISS(  1)  .EMISS(NLIHC)) 

011380 

BUFFER 

OUT(MFILE,l)  (FF(U  ,FF(NllMO)> 

311400 

BUFFER 

IN  (LFILE  ,1)  (Xin(l),XI0(2)> 

0117B0 

BUFFER 

IN (KF ILE • 1) ( XIO ( 1)  , NL  A YER) 

011800 

BUFFER 

CUT (MFILE, 1) f-IO(l) ,NL*TER) 

011020 

5 35 

BUFFER 

IN(LFILE,lHOvlP,NLIHO) 

011950 

BUFFER 

IN (LFILE, t) (OLOEM(l) ,OLOEM(NLImO)> 

011970 

BUFFER 

IN (LFILE, 1) (OLOTR< 1 ) , OLOT R( NL IHO) ) 

..1199C 

BUFFER 

IN(LFILE.l) (OVlP,NLIMO) 

012423 

BUFFER 

IN(LFILE.I)  (0L0EC.(1),0'.QEH(I4LT"0)I 

312480 

BUFFER 

IU(LFILE,1)  (OLCTP(l)  ,OL0TR(NLI*‘O)) 

012530 

BUFFER 

IN(LFILE,1)  <OVlF,NLIHO) 

.13.6J 

BUFFER 

IN ( LFILE,  1) (0LCEN(4) ,OL OEM(NL IMO> ) 

01 31  JO 

130 

BUFFER 

IN (LFILE ,1) (OLOTR(4) , OLOT R (NL IHO) ) 

013130 

150 

BUFFER 

OUT (MFIL  E, 1)  ( V IP » NL  IM) 

313420 

BUFFER 

Ol-KMFIL1;,  1)  (NEWEM(l)  ,NEWFM  (NLIMI ) 

*1  34*,  j 

BUFFER 

CUT (MFILE , t ) (MEWTR(l) , NEHT» (NLIM) > 

013460 

360 

TUFFFR 

IN(KFILE,1)(V1P,NLIN) 

Cl  36  33 

BUFFER 

IN  (KFILE,  1)  (TP(1)  , TR  (*IL  IN)  ) 

.1 3653 

bur  re  r 

IN(LRILF,1>  (Xin(l) ,XIO(2l ) 

014570 

BUFRC9 

IN(KrILt,l) (XIO(t) ,NL*V£») 

014590 

SURFER 

OUT (MFILE ,1) (X IO ( 1 > ,NL»YER) 

, 1<*61  j 

535 

BUFFER 

IN(LFILE, 1) <0V1P,NLIM0) 

C14740 

BUFFER 

IN (LF IL  E , 1) (OLOEK(t) ,0L9rM (NL IMO) ) 

Jl 4763 

BU'FER 

IN (LFILE, 1 ) (OL  OTR  < 1 ) , PL  BTR (NL IM0) ) 

014700 

9UF  R£R 

IN(LFILE,1) (0V1P,NLIN0> 

Cl  52  20 

403 

BUFFER 

IN  (LF IL  E , t ) (OLOEM (l) ,0LOEH(NLlNO)) 

0152S3 

BUFFER 

IN (LFILE, 1) (Ot  OTR( 1) ,OLOTR (NL IMO>> 

91  53 30 

BUFFER 

INUFILS,  1) (nVIP.NLIMO) 

015803 

buffer 

IN (L* ILE,  1) <0LPEM<  4),0L0EM( NLIMO* ) 

315040 

42 

BUFFER 

IN (LFILE, 1) ( OLCTP (4 ) ,OL OTR ( NL IMO> > 

01*060 

150 

BUFFER 

0UT(MFILE  ,1)  (VIP,  NUN) 

01 ol 5u 

BUFFER 

OU T< MFILE, 1) (NEWEN(l) ,NEW£M (Nil*) > 

31617C 

BUFFER 

OUT (MFILC ,1) <NE NT R (1 ) , NO Wf R ( NL IN) ) 

C 16194 

360 

BUFFER 

IN(KFILE,1) ( V1P,NLIM) 

016390 

BUFFER 

TN(KFILE.l)  (TR(1),TR(NLTM)» 

11641k 
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TABLE  B.2.  SAMPLE  RUN  DECK  FOR  TEST  PROBLEM 


74/74,  OPT-l 

FT N 4.6*425 

3 3/ 14/7  S 16.43.20 

CLOUGH, CH577fiO,T48. 

FASCODEL 

Z <, 

PLOUGH 

ATTACH, TAPE?, CL0U199 3TO2210 , IO-CLOUGH,M»»l . 

REQUEST, TAPE12,»PF, 

A TTACH,  HIRACl , ID-CLOUGH, HR» 1. 

FTN,I»HIPACl,Sl,8»lGO. 

ATT  ACH,H  IRACD,  TO* Cl CUGH, hp»  1 . 

FTN:l«HIRACn,Sl,B«Of.n. 

ATTACH, VGO»FA  SC  ODE, t 0*CIOUGH,H\ »1 . 

PEHINO, OUTPUT, 

HAP , PART • 

LOSET, PRESET- INOEF. 

LCAC, VCO . 

LOAD,  ISO. 

IOA  C,  OGO  • 

SEGIOAO. 

EXECUTE, FASCOOE. 

E09 

FAST  TPEE  FASCOOE-1MIT ACV.HTRACL.HXRACO,  APS MRG , E MINI T , E MUF ,£ MPOWN) 
FASCOnE  GLOBAL  HAIM, HEM 

HIP AC V INCLUDE  SM APE l, SHADED, VOT CON, BPFI LE,C ONVF  HW  .PANEL ,HOL  EC  .OVRFAC 
HIP ACL  INOLUOE  SMAPEL,  ROFILE.CONVFNL, PANEL, HOLEC .OVRFAC 
MIPACO  TNCLUOE  EHAOE0,RCFILE,C0NVFNC,PAN£Lr**M0LEC,0VPFAC 
END  PA  SC  OPE 

CO® 

1 APPPflX  TO  KYLE  PEPOPT...  . HIOLATlTUOt  SUHMS9 


1 L 

1. 

0 0 

1 

2 395. 

2105 . 

273. 

10 

525.0566 

264.3769 

4. CO  KM  TO 

x.OC  KH 

9. 38flF*Zl 

1 . 43  PE  *21 

2.611E*17 

1 • 195 E*tP  3.201EM7 

6.931F*19 

8.942F+23 

3.334F+24 

332.0X41 

242.7476 

7.01  KH  TO 

11.30  KH 

1.649E*21 

1.284E*21 

4. 987EM7 

1.3S9EU0  ?•  91 /r* 1 7 

a.’24F*13 

8.149F*23 

3.03/1*24 

197.  04  44 

222.0735 

11.00  KM  TO 

14.00  KH 

0. 81 5E*19 

6.262E*20 

9 .351EM7 

5.313E*17  t.423E*17 

3 . 0 35  F ♦ 1 J 

3. 97EE ♦ 23 

1, 402E*2h 

124,0100 

216.CC00 

14,0"  KH  TC 

17. 3C  KH 

7.6l2E*t5 

4 . 04  3E*2  0 

7.052E+17 

3.429E*17  E.109EM6 

1 . 960 E*  1 9 

2 .566 1*73 

9.5632*23 

77.2771 

216.9193 

17,03  KH  TC 

2C.r0  KH 

5. 045E*18 

2.509E*20 

1.0  9l£*l9 

2.129EU7  5 « 7 "2E ♦ 16 

t .216c*13 

1 ,59JE*?3 

%936f*  23 

36.4133 

223.8«96 

2C."0  K.-  To 

30. 0(3  KH 

1, 349c*19 

3. 204E  *20 

x,33 1E*13 

2.786E*17  7.462t*16 

1.592C31? 

2.085F*?3 

7.  7 7 :f  *2,3 

0.3432 

242.9639 

30, OC  KH  T C 

4". C C KH 

4.904C1A 

7.053£*19 

1. 254C13 

5.902EM.6  1,6332*16 

5 . 4 19E  *1 ’ 

4,‘.77F*22 

l.b  6 )£*23 

2.44904 

26-».2r7n 

4“ , C 3 KH  TO 

45, 30  KH 

4. 913c* l 7 

1.110E+19 

1.533EM7 

9,484E*15  2.541£*15 

5. 42.E+16 

7,  C9 7E ♦ 21 

2.649E.-J 

1.1554 

267.2521 

45, CO  KH  TO 

54, Of  KH 

2 .59  CF* 17 

0 , 43  7E *19 

4 .997E»ift 

7.19eEUE  1.917E*  15 

4.  0 9CF*U 

5,355  E*21 

1 .996  E *22 

.2742 

255.4215 

94,00  KH  TO 

1 1 J , 3 0 KH 

4.745E*16 

4.050E  10 

1.174E  16 

3.464E  ,5  9 , 276 E 14 

1.979E  16 

2.591E  21 

9.64  9E  21 
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V-w' 


XID 

IHIRAC 

I EMI  T 

I LOOK 

I PLOT 


TABLE  B.3.  INPUT  DATA  DESCRIPTION 

7 element  array  contalnint  text.  The  text  is 
printed  with  FASCODE  output  and  is  used  as  a 
header  for  execution  time  plots. 

1,  Voigt  approximation  to  line  shape 
P,  Lorentz  approximation  to  line  shape 
3,  Doppler  approximation  to  line  shape 

0,  no  emission  calculation 

1,  emission  calculation  performed 

0,  emission  calculation  from  ground  to  space 

1,  emission  calculation  from  space  to  ground 

0,  no  plot 

1,  on-line  plot 

2,  off-line  plot 

3,  CRT  plot 

4,  graphics  plot 


ISCAN  Indicates  parameter  to  be  plotted  (use  if  IrLOTfO) 

1,  absorption  coefficient  YY=YY 

2,  log  of  absorption  coefficient  YY=AL0G10(YY) 

3,  transmission  YY=exp(-YY) 

SECANT  secant  of  the  angle  between^he  1 ine-of-sight 

and  the  zenith 

\ 

VI  (lower  limit)  specify  the  wave  number  range  for  the  output 
V2  (upper  limit)  ^ * 
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TABLE  B.3.  INPUT  DATA  DESCRIPTION  (continued) 

temperature  (°K)  used  to  compute  a boundary. 

If  TB0UND*O,  no  boundary  Is  computed  (currently 
used  for  looking  down) 

number  layers  + 1 


Pft.VE  average  pressure  (mb)  for  the  layer 

TAVE  average  temperature  (°K)  for  the  layer 


WK  array  containing  the  molecular  column  densities 

p 

for  the  layer  (molecules/cm  ) 

h2o,  co2,  o3,  n2o,  CO,  ch4,  o2 

**★***★★* a***************************************************************** 

IF  I PL0T=1 


PLOTID  3 element  array  containing  30  characters  of  text 

used  as  a beginning  and  ending  banner  on  pen  plots 

VI  th«  range  of  wave  numbers  over  which  plot  will 

^ made  (cm"*) 

XSIZE  size  of  X axis  (In) 

DX  number  of  units  of  X per  Inch 

YMIN  range  of  Y-axis 

YMAX 


wwsc 


the  constant  k in  the  expression  y=e 


TABLE  B.4.  LOAD  MAP 
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TABLE  B.5.  TEST  CASE  OUTPUT  (continued) 
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TABLE  B 5.  TEST  CASE  OUTPUT  (continued) 
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TABLE  B.5.  TEST  CASE  OUTPUT  ( continued ) 
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TABLE  B.S.  TEST  CASE  OUTPUT  (continued) 
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THIS  PAQt  IS  BBST  ^U^1^f|PlUCTl^M* 
TOO*  non  FlIKWiiiHii)  10  WQ  — 


»rt 

0 


M «v * 

4 » 
hM 
4o«t 
0 m * 


.,  mtfctwu>wuuib  ui  i 

0 W *1  v»  w0  *>  -*  4 

m00,#K»vOvP*4«J* 
n^»<*n4fstu>us<j  j 
V ^ ^000000 


«K.0tt*4*W0O'  0 

»•  «* 


Sl>  M UN  «k  0 • 

f • 0 

*0  (T 


«a0«0OWGI0'~ 


*» 

J 

ftl 


UJ 

0 c 

« c 


• I I I * I ••••••  • 

IbUU  IM1'  I4<  U.1  Ui  U Ui  IU  41 

■-*•»»«*  0 «NNn4\f«|fi  0 


K)t  K 00« 

4«4*«nrtNNNN(4 


*.  ft-  K 

4 5h 


O'  * (ft  «•  0 
M r~  * 

e»  < 

O'  0 I 


«>  *»  <D 
0 0*  0* 


m00®000000000 

sK««0V^f4NN.Y«  J 
rr^ir^raoiiV'A*  * 

RRRSERSSSrSSS 


ts 

Ui  O 

xte 

SS 

Ul» 

«ts 


W I 

* 


* 

O' 


» 

i 


iii 

Ch 

• o 

S* 

^ V 

ss 

• k. 

• UI 

• 0 


> o 

ft-  O Ui  UN 
ft*  I* 
► ^ 3 0 

< 0 I • 

j«Cn 


c: 

Ze 


C 

9 


4 

** 

V 


s - s * 

8*4  « 0 

• 0 


V 

V 

t-  o 

u* 

o* 

*g 


r 

c 

X 

IT 


1U 

*-  w 
*•  « 
9 

-4 

e 

> 

i 

i 


CN 

*<  -H 

s* 


♦♦♦♦♦♦♦ 


o 

* 

a 

H 

u. 


& S u*  i i i i 

• -i  bkk  if  n 

K © 0»«  0 IT  •* 

« f 0 0 0 *4  0 

P4  * «•*  • • • • • 


& s 


9 

O 


u 0 

0 0 


a »i 
M 0 

0 «• 
►*  • M 

* 8 r 

u w r 

W ft  Ui 

0 ft  h- 


* - 

S CN4CC  4M  ~ 

* s? s°aus°  ? 

5 5 

* > 

u o 


0 

0 


* 

£ 


107 


S^WSSaraSS 


40U«nvS(C*-ll  « .147? 


TABLE  B.6.  TEST  CASE  OUTPUT  (continued) 
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TABLE  B.5.  TEST  CASE  OUTPUT  (continued) 


m 


THIS  {“AOS  iS  BEST  QUALITY 

nm  ccn  fUKtushiOJ  xottifl - 


•> 


N 

IA 

UJ 

4 

I 

C 

o 

* 


W1 

a 


.1 « 
UJ  fu 
* <vj 


V K 

o «. 

M to 
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* 

t 

o 

u 


OK 
« *4 

at 


* 

♦ 
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»*  * 

s : 


nc'^ermo«^<»«nn,er  3 

I ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

tw  y id  w ui  unu  iv  ui  l ui  vu  w u u< 

# n tr  ^ a a m k 2><n2u(r«Q 

<AV'^^CHNKfniM^MK^e»4 
U'K  ■ J A j 4 *4K  ,**4  *■*.>*>*»*'•  -» 
^v4^ir&»*4fvj<y«ia9'«£M>K«>«>ru 
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er*  ■* 

«•  *4 
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USAGE  NOTES 


1.  The  user  must  structure  the  input  atmosphere  such  that  the  ra  io  of  the 
previous  DV  to  the  present  DV  is  acceptable  to  FASCODE,  (see  the  sample 
input  data). 

2.  The  program  is  written  for  execution  on  a CDC  6600  computer  system  at 
Hanscom  AFB.  Accordingly,  the  control  cards  Input/ output  statements 
and  plot  options  need  to  be  tailored  to  the  user's  installation. 


THIS  PAGE  IS  BEST  QUALITY  PRACTLCAWi 
mm  COPY  FURSUSHED  10  DDC - 


PROGRAM  FAS**0flE  74/74  0PT«1 


FTN  4.9*4Z«  03/14/7* 


1 


9 


19 


19 


ZD 


?9 


30 


39 


40 


49 


90 


99 


90 


6? 


70 


PROGRAM  -'r*SCOnr»lHPUT«133/'lJ05,OUT7'  *130/1300, TAPE94INPUT 
r .TAREP-OUTRUT.TARE?,  TARC10,TAREli  'ARE1ET 

r 

C EAST  SIGNATU°E  COTE 
C 

CCMNON  CC*?TRtS»0,S* 

COMMON  77<3900»  .srifoo»,vsrtzzsi 
COMMON/MAIN/  KFIL*,FP»NEI.,*fl,TEMPO 

COMNON/NEN/  NM0L,XIIM7»  , SECANT  »P  AVE.TAVE  fPHOLIO  ( 7>  ,MK  t 7)  , 

C OV.Vl, VZ.NLATER 

COMMON  /8UFH0/  FILMO < 1C > ,9HOLIO<10>  , HOllNO (1 0 1 > VLINLO*  VLINMI 
C 

NMOL* 7 
ITYPE—1 
7041013  • 

Tfnfo*ZB6. 

XRANEL*0 
LINFIL*  3 
HFUf4lZ 
LFTLE*11 
KFILr*l 0 
NLA7EA4P 
p»INT  *00 

REAP  910,  (XIO(I)  ,141,7) 

C 

C TMIFAC  4 1 CALL  MTRACV 

C TMIFAC  « ? CALL  MIFACL 

C IMIRAC  ■ 3 CALL  HIPACP 

C 

C IF  IENIT  .EC.  1 PROGRAM  MILL  COMPUTE  5MISSIO-. 

C 

C ILOOK"  1 L COKING  FUM  SPACE  TO  GROUNO  FOR  EMI'S  TCM 
C I LOOK*  0 LOOKING  rRO"  GFCUMO  TO  SPACE  FOR  EMISSION 
C 

"EAP  9Z0,  IMIRAC , IEMI T,  ILOOK, I SCAN,  TPLOT 
IF  IIEMIT •EN.O)  XIO  <7)«10M  A»S  CP'FF 

IF  < (I'MIT.EO.l  > « AMO.  JILnOK.EO.OT  > XT0I7I*1IM  EMISS  UP 
IF  < (TEMIT.E'M  ) .ANO.  < ILOPK.EO.  1 > > XIO«7)*IOM  ft*ISS  OM*l 
P-*IUT  9 30  i XIO 

FPINT  «0,  IMIRAC,  IFFIT,lLOOX|I  SCAN.IPLOT 
RFAP  940, SECANT 
P“INT  990. SECANT 
PRINT*,  "IMIRAC  « ".IMIRAC 
C 

C T"OUIIO  IS  A TEMRE4ATURE  MMICH  MILL  PS  USED  TO  COMPUTE  A 
C FOUNDRY.  MMEM  T9OUN0  ■ 0,  NO  BQUNPRT  MILL  01  U5EC, 

C 

• FAP  940*  V1*VZ 

RRIMT  »,  " VI  , WZ  “ , Vl.VZ 

RE  AT  «4p,  TROUNO 

PRINT  •,  " T BOUND  " , T90UN0 

C 

REAP  9?a,MNlAYP 

BUFFER  INUINFIL.i)  (FILMHI1), VLIMMT ) 

IF  I UNIT  (L  INF  ID  ,EQ* tt  STOP 
PRINT  930,  FUMS 
PRINT  009 

PRINT  090  , (BMOLIP<l»,MOLINC«I».l4l*NPCU 
P*INT  070 , VLINLO, VLIMHI 
IF <Z» (N4LATR/Z) ,NE.  NNL AYR)  GO  To  110 
NSTORoMFILf 
MFILF«LFILE 
LFIL»4MSTOR 
110  CONTINUE 

I'  UEMIT.EO.  P>  XFIIE-MFTLE 

00  ZOO  K*1,MNLAVR 

POINT  9 Of 

IS  PALI* 0 

REMIND  LIMFIL 

BUFFER  IN  «LlNFI*.lt  1FILHO  (1  > , VL  INHI) 

IF  (UNimiNFIt)  .£C.  0 > STOP 
MLAYSRoML  AYEP*) 

PRINT  *,  “ MLAVEP  ~ , MLAYC' 

IF  (K.EC.NNLATR>  XIC(9J4 lOHLA'T  L»TER 


000101 
00011C 
OliOIZO 
Ot  Clio 
000140 
C00150 
000190 
000170 
OCClAi 
000190 

ooozoo 

otozio 

ocozzo 

000Z3D 

0C0Z40 

GC0Z90 

009Z90 

ccozro 

OCCZSu 
or  oz90 
000300 
000310 
0C03Z0 
000330 

Pt  0340 

SL  1390 
00  0390 
OtilTC 
GCG30J 
Ot  0390 
000400 
Cf04ir 
0004Z3 
000430 
Cl  C44C 

r"045c 

ori.46C 

C0C47C 

'to-.sc 

ri  .-9 

000900 

:( cot" 
:c  csz^ 
000930 

;i  C94-, 

if C990 
000990 
:i  C97C 
CC09A0 
JCG99C 
ci u9o: 
CL (910 
or  zb?: 
CLC9T0 
PC  0940 
PC  095C 
000990 

eroirc 

CC  0990 
01 0990 
OCOTiO 
000710 

n 07zo 

0 C 07  30 
000740 
OC  C79C 
oro79c 
000770 
000700 
OC  0790 
OGuOOO 
OSORIO 
OCOOZC 
CPOISC 
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SHIS  PA(ZE  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPY  FUKMSHJSD  IX)  DOC 


MMMII  wrtm  T%n%  OFTM 


ftn  m*mi 


ft 


•I 


tt 


tt 


M 


111 


ill 


111 


ttt 


III 


ill 


tit 


tit 


tM 


l.t 


wto  tM,nitiim 

fttAP  tM.INKnKI.XKn.Tl 

in  oiooccv 

ot  me  o»«m  vt /•  t • ton  t i ttwt  t r t •««) 

UIU*ll,tll*NII«tOI 

IIUMKIi. 

iDin<i<iitMi«ii'ii.iam«vniMoir  iTtvt/avNtni 
Mtt»M.t*(  nut  * i out  ( tjitmut  * i.'toun'toin  > » 

XF4  tHXMC«(Q. tt  OVMVttt/.. 

XMlNTRtC.n.tl  OHHHI/i, 

III MlHtO.tO.lt  OnMHI/l. 

It.  COtTTKUt 

01  • ",0V 

It  INitVtK.Nt.lt  00  TC  tit 
C ««  It  IttUNCO  TO  tC  .IT.  t 

mu*uotuiovi  - ?. 

Mu*tt.»*ntti 

tt«i(M/teiu  ♦ t.t 
ov*tcw.*no«Titov) 

to  TO  ill 

tit  TVFt«OtOM/«V 

ITITTtt.lt. t.t  to  Tt  tit 
Ttn«l.  l/TVOf 

tmu«t 

in  FKlNT*,«COI#UTtO  OV  Ottott  NOOIFXCtTTON  ■ *»0T 

0 

itmtr  .tt.  t.tt  to  to  ttt 
ItITTtC  .or.  t.tttl  to  TO  tit 
ovaotcot 

TTTt«M 

tC  TO  ttt 

tit  ITTtf«t •/ |TVFX“1 « tl *t  «t 

xrutMtu  .to.  «)  ot»aooy»ti.otT«xTTtn/ttotT«tTTtt*n 
irasntu  .to.  i»  OT»otoi»*tt.otTaTTH*ti/no»T»xT»tti 
to  To  t»t 

ttt  CONTINHt 
••XNT  ttt 
STOt 

t»l  F*INT  •.OT.TTTtf 

c 


c 


ttt 

ttt 

c 


0 

SOI 

0 

Ilf 


lit 


ttt 

lit 

r 


IMlMXMC.n.t)  GUI  NIMCV 
tMINXttC.ro. II  ctlt  NlttCL 
XMIHXttC.CO.TI  C*U  HXMCO 

NtTSvt 

xr  ixfNn.w.ti  to  to  ttt 
IF  lHttm.tT.il  to  TO  ttt 
to  TO  Itt 

cut  winno  iiTTNttXtFttt.xeioNt.tme.NFitn 
to  TO  ttt 

CONTINHt 

Xt«XT»t  TO  CMC*  THIS  STtTCNINT 
XF  IItNHt.NC.ll  STOF  " ISNUt  .Nf.  t “ 

XF  <NltVtK.tT.il  10  TO  190 
CUt  CNXNITIXtOOK.TtOUNO.NFIttl 
CO  TO  Itt 

IF  lltOOt. CO. tl  CUt  CNUF  INFTI,tTTFt,tFXtl,HFIt*l 
XF  atOOK.CO.il  cut  INOONN  <NFTS,ITVFC,SFH.l,NFIlCl 

CONVXNWF 

HSTORHNFItC 

NFXttHFXtt 

IFXIC«NST0N 

NCNINO  tl 

Mint  tt 

xoiOFi«Kotnct 

KFINCIM 

KFXtC«tt 

CONTXNUl 

OIKTNt  It 


C 

ttt  tF  UFtCT.NC.  t I CUt  TFtOTlTFtOT.NFUII 
lit  STOF 


tlFtt/lt 

ttttM 

line 

ccmo 

ttttM 

KlttS 

ton  to 

0 DSttt 

tottto 

tttttl 
tilt  It 
tttttl 
tttttl 
tttttl 
tottrt 

tottto 

ttotto 

tottto 

ItlJil 

0(1410 

te  to  so 
ttts.o 

tottto 

tettto 

oeitro 

0C10I0 
CtlOtC 
OC 1100 
OCtllO 
0C11I0 
otitic 
tlllll 

tottto 

trtttc 

00 1 1 TO 
C 0 1 1 In 
crutc 
ociroe 

OOlJK 
oriJto 
0»  tilt 
tl II. c 

acute 

ootm 

Otitre 

atiltc 

totitt 

ttnoo 

oitna 

ootm 

OllllO 

ocis.c 

ormo 

tint' 

tcisre 
toistt 
oemo 
tr t.oo 
ott.to 
oci.ro 
cci.io 
tot. .8 
ODl.tO 

oct.ts 

cst.ro 

01  t.tc 
toi.to 

totttt 

osmo 

ocitro 

ttiiii 

Otlt.C 

octtto 

tctttt 

oettTC 

tri.it 
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JHIS  PAQS  IS  BEST  QUALITY  PTU.CTICAW4 
I^OM  OOfX  JUHJUSHBD  XO  DD.G  ...  — 


rognip  mwom  o*t»i 


ftn  4.4»iu 


DMWM 


Ml  F0NN4IIIMU 

««  roixtniHSi 

no  fonhitiiiuii 

its  rOKMItlMII 
MB  FO4NAT<lH4,llAt0> 

mi  foxnituofii.t) 
mo  fo»h»t(»o  sec»HT  • •,*!*,»  » 

lit  FC4M4TI4X .41*  3N  > .11) 

(F6  FCI»N*T(*I  LONIST  LtNf  • •Ftl,S,»  HXIHMT  LI  HI  • » FIO.S) 

410  F0KM4TI*  THf  NIT to  OF  010  OV  TO  NfN  OV  (XCtfOS  t.O  • ) 


001MI 
001100 
111410 
001410 
1(100 
00 1440 
001440 
01 100 
0014  ''0 
001440 
OC 1440 

loirao 

iiifio 

OCIFfO 
SO  1 F 39 


SUN  ROUT INC  MI *»rv 


F 4/F*i  C»T«t 


FT  M 4 . *4  4(1 


OS/14/FO 


SUIROLTXNC  MIMCV 


CILCUUTfS  MONOC  HNOMATT  C I4S0NMT  ION  COfFFIClOT  FOI  SI4GU  L4VIR 


U4C4  4FF40XIN4TC  VOIGT  OLGORITMM 


COMMON  GNU  1(401 . S((40l , *LF*0  1(40)  ,(M(fS|>  .NOL(MC)  , 

C tFOFTMtdll  ,N(C»lM(IS0>  ,HT»1 1(40)  , I(C T»((*0  I 

COMMON  FF (3400)  • SF  (40  0)  . VJF  1 ( ?4) 

COMMON/M4IN/  XFILf  , XMMFL.FO  .TFNNO 

COMMON /NCN/  NMQt , XIO (F) , 5KC4NT. POVf , T*  VO  , MNOLIDI F) ,KX(F>  , 

C OV. VI . V( .ML 4TIN 

COMMON  /XSU4/LIMIN,  ILC.  XNt . MOTvVTOF  . VFT  . I (OF  .IFONIL.ISTOF.IOATI 
COMMON  /MXX/  NNF.OXF. HF.NNS.OXS.NS.NNVS.OXVS  »NVS .NFNAX .NSN4X .NV4M 
14X 

COMMON  /SUI1  / M4XF.MXS.M4XV4.NLIMF,  NLIMS.NL  IMV1.NLO.NMI  ,OVS  ,0VV5 
COMMON  /XNINCL/  VIF.VFM.OVF.NLIM.NIHTFT.NFTS 
COMMON/  XT  INC /TXMf  ,T  IMDOF . T IMC.NV,  TIMFNL 
COMMON/ VOICOM/»VN*T((01)  .CLO(IOl) 

OtMCNSICN  »LFCO»(T)  ,»LF01(F),g(F)  .SCONIF) 

OIMfNSION  XF ((41 1 , XI (111 ) . X VS  1(41 ).X0((41> 

06T4  NNF  / 4/.  OXF  / l.lt/.  NF  /(II/.  NFN4X  /(It/ 

0*T»  NWS  /14/t  0X4  / 4.00/.  NS  /(I!/.  NJN4X  /(II/ 

0»T»  NNVS/44/.  OX VS/  O.lt/,  NVS/II1/.  NVSN4X/I41/ 

04T4  If NTtN/0/,  LtMIN/140/,  NSHIFT/St/,  NLlMF/1401/,  NMTS/  0/ 

04  T*  SUIIO  /ION  VQTCT  / 

IF  (TFNTT4.Nf.il  GO  TO  10 
IIMT14»1 

NLINSm(NLIMF/4I*1 
NLINVS*  (NLIMS/4)  41 

NOT!  (OXVS/OXF)  IS  14  4N0  (OXVS/DXS)  IS  4 
N40UN0*  FLOAT  (NNVS1  *(0  XVF./'JXF) 

NAXFRK«  INFfNSOUNO 
M*XSm  (M4XF/4I*! 

NIXVS«(M4X3/4!41 

OKI  VO  ICON 

C4LL  •MMKXF.XS.XVS) 

C4LL  (NA4f0(X0> 

C4LL  NOLtd  1, MMOL  10, NMOL, TfMFO,T«VI,M»,F»Vf, SC04,»LFCC*.4LFD1I 

MINT  444 


001F4I 
001F40 
001F4I 
001FF0 
QC  IF  40 
001F40 
OC  10  00 
Ot  1410 
0014(0 
001630 
OC  1440 
0(1*40 
001440 
Oil  4F0 
001440 
9114  90 
001400 
101110 
0011(0 
001130 
001440 

oeiiio 

00 1440 
0I14F0 
001410 
0C144O 
00(000 
001010 
eoioio 
00(030 
00(040 
10(040 
00(141 
•4C4FI 
41(444 
•4(444 
04(101 
14(110 
00(1(0 
00(130 
10(140 
001140 
0*1140 
OIIAFI 
OOXIOC 
00(144 
mill 
•4(414 
MON 

•inti 

••(CM 

•him 
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THIS  PAGE  IS  BEST  QUALITY  FRACTICAKL1 
FROM  COPY  FURNISHED  TO  DD.C 


MMNwrni  kimcv  fwn  cpt.i 


ptn  ii/u/ti 


ft 

if 

if 

» 

rs 

•a 

tf 

90 

If 

no 

Iff 

Ilf 

US 

Iff 

its 


MINT  STf,  SUB  10 
T INPOPaTINCNVaT INPNLai. 

MINT  994,  <KXO<X)tX«l(r> 

Rtf  SBC  rfILE 

I fora  I 

MINT  fit,  «rxu 
MINT  Oil,  SICINT 

mint  on,  Mve.Tift 

MINT  Oft,  IMtiniW,W<IH,Ml,mtU 

MINT  94t,  90,01,01 

000*90 

OOlalOXVOm  *00 
DV«tafom/BKr)«00 
tOUNOrOiOITIMOUNOI  roo/f. 

•OUNOOafOUNO/t. 

MINT  OOt,  OOUNOt 
AL(rMAt»90UN0/PL0AT<NNVlt 
ML 0*N OOUNOT' . 

NHIaNLINPaNSHlPT-i 
00  St  lat,*»xr 

m rrutat. 

DO  It  lal.MXS 
It  SPllfM. 

00  7f  jal,MXOS 

rt  osrufat. 

c 

$ UNITE  fKPlLEt  IXlOU)  «lal,7t, SECANT, PAVE, TAVC,(NN0LI0tm  , Hal, 7t  , 
t * <NMI»,Ial,NNOt>, 00,01,0? 

supper  outikpile.I)  ixidiii.nlavcri 

ir  lUNIHKPILCI.fO.0t  STOP  -99- 

c 

OrT»Vl-?»*tOUNO 

OOOTaOl-OOUNO 
OT  ora  or*  NOUN  0 

c 

XKTOat,  ISS1*TENP0 

XKT  ■ 1,1091* TINE 

VKTFtO  (1,/XKT)  - C1./XKT01 

CALL  HOLEC 12 .HI1011D, NMOL  .TEMPO  , TAOE*  NO,  PAVE,  SC  OR  .ALFCOR , »LFO  1 i 
00  80  Ml, NMOL 

•0  LMN1*NK  CMI*SC0R<M1*SECANT 
ICATaO 

SUMALP«0. 

SUNTETaf, 

NCHNCaf 

C 

OO  CONTINUE 
C 

100  CONTINUE 
C 

CALL  SECOND  ITTME01 

if  ucDr.Nc.ot  co  to  iso 

CALL  NBPILStCNU»S«ALf A8,CPP,nOL> 

CALL  SECOND  I TIME) 

TINNOFaTINNOMTTNE-TINEO 

C 

IFtlEOF.ItE.II  CO  TO  lit 
C 

C HOOiro  LINE  oata  ro»  TEMPERATURE, PRESSURE,  *nc  COLUMN  density 

e 

oo  no  i*ilo ,ihi 

MaMOMI) 


1(1211 

082270 
082280 
0722  SO 


002320 
Of  2330 
082311 
0823S0 
082310 
882371 
012310 
08 ESSO 
082100 
002118 
082120 
882430 
88  2110 
802490 
072110 
882470 
082180 
082490 
CI2S00 
882S10 
0C2S20 
002S30 
0C254O 
0C2SS0 
DC2SI0 
08  2S70 
l "2580 
hi  2590 

or  2iof 

082613 
08  2120 
otecso 
002880 
012650 
08  2668 
002670 
002680 
00269C 
DC273C 
002710 
08  2720 
OC  2730 
0C2740 
002/50 
082760 
002770 
002780 
OC  2790 
002800 
802810 
002820 
302830 
002840 
002850 
002860 
082870 


CFCPTP(l)aStIt*UtMl  002880 

ir  ceropTMdi.LE.o.)  co  to  no  ocze90 

ICNTa ICNT41  002900 

ALELaALPAO  (I)PALPCONtNl  002910 

ALPADaSNUlIf  *ALP01(P)  082920 
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rxMMINaCXM4-lTL.il 

COll  SfPONOITIMCI 

MOINT  *,  - INUM,  MM  lit*  LMILt.  KMItt  MM1L  t,  LMItr.KMU  t 

MMIMT  220,  Tint 

im«t 

•MINT  THIS  OUTMUT  IS  MMITTCN  ON  Mitt  MMILt 
C 

KNlaXMONCL 
tO  MtMINC  MMILt 
NtMlNO  LMILt 
NtNINO  MMILt 

NUMMfM  IM  ILMKt.lt  ixiniit.xinitn 
IM  IUMIT4 LMItt)  .fO.  0 I OTOM 
OUMMCM  1M4  KMILC.1I  IXIDtlt.NLOTtNI 
IM  4 VNtT  4KMIL  f t .tO.  ft  ST9M 
•OMMtM  OOT4MMILM.il 4X10411 .NLITfftl 
IM  II'NIT  IMMItf  I .CO.  ft  OTOM 

oo  so  x*i.o 

CMOtVtlKtaf  . 

Of  TM5»Vt4K)»0. 

»TTPt«ITVPt 

1M4ITTMI  ,Nt.  01  GO  TO  0S9 
C 

C 1/1  NOT  10  ONLY 
C 

JUGaS 

JOC«f 

CO  TO  900 

ISM  NUMMMN  IN4LMlLt.lt  I0V1M.NLTM0) 

IM4UNIT  4LM1LCI  .CO.  0)  STOM-f- 

04IMMCM  IN  4 LMILt  .l)IOLOCM(l)  .OLOEMtNLIMOt  t 

I*  4UNIT  4LMtL  ft  .10.  0)  JTOM-S- 

BUMMf  * lN4LMILt.lt  4 0LCT4M1I  .OLOTNINLIMOt  t 

IM44INIT  ILMILMt  .CQ.  Ot  STON"*“ 

00  95  0 la l.NLIM 

NCNCN 41 |a OLOt M 4 1 1 ♦ 4 1 . 0-T94 1 1 1 •■••0L0T9 4 It 

NCMTN4IlaTN4It*OLMTN4Il 

TtJTf4t*TtSTf4toOV* 


onooo 

01 10  90 
011000 
011010 
Kioto 

011030 

0119*0 

011300 

011000 

ousts 

011000 

ousso 

onooo 

onoio 

onoto 

01 10  so 

01104.0 

onooo 

onooo 

onoro 

011604 

onooo 

onxoo 

Gimp 

011X20 

onxsc 

011X09 
011X01 
011X09 
OllXXt 
01 IX  09 
011X00 

onooo 
011010 
011020 
tnosc 
0110*0 
’Ilf  0.1 
011**4 

111  ox; 
CUTS. 
01109? 
Jim: 

ohmic 

111924 
01 19  SO 

oimc 

011990 
011900 
0119X0 
oi  m? 
9119941 
012000 
012119 
C12020 
oitosr 
0120*0 


IMTM5TC*  ,LT.  CNONCC)  <<0  TO  ISO  D12j9C 

C0MCNaCM»N6t-M00  912.00 

•|aDO0ItS4LTtSTI- 4000ItS4LTf  STI -9001  CS4LTfST  Mitt  *COHON  0120X0 

CM«NCfaCH»Nat*0V9MM  912000 

IM 4CNONCC  .Gt.  0004  1.0)  CO  TO  900  412090 

PC  TO  SSI  012100 

900  LTtSTalTCSTM  012110 

I00a|004l.l  012120 

SSI  COttTINUt  1121  SI 

•0  TO  1U  0121*0 

ISO  IKNlalKNLM  912190 

IM4IKW  .tt.  KMLt  r.0  TO  SCO  012100 

CO  TO  210  OtllTO 

C 0121044 

C Ott  MTIOS  fXCfMT  1/1  012110 

C 012200 

*St  ttalTTMO.l  012211 
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*ffis  PASS  rs  BJST  QUALITY  PRACTISAfclii 

ZRQU  QCVPX  PURmsHED  IOLDO 


subroutine  ehuf 


74/ T4  OFT-1 


FTN  4.6*428  OVU/TI 


i 


I 


f 

i 

[• 

! 


t 


TT 

•I 

II 

H 

n 

in 

in 

ui 

in 

120 

its 

nt 

131 

1M 

III 


wsnfnnronn 

JCC«1 

oo  ii  jnaiiirm 
AFI-JFG 

MINT  -i“  F-  -,F 

C 

C ••••*  THE  FOLLON  ARC  TNI  CONST  MTS  FOR  THE  LAGRANGE  I FCINT 
C INTCRFOLATXON 

e 

IKJFSt  — F*«F-1.  II  • CF-E.  |)  /I.  I 
UIJKI*(H«M>I)*IP*M)*I.I 
l3UKI««l*(M.lt*IM.n*l<l 
•l(JFCI*F*|M*t“l»ll/l«l 
II  CONTINUE 

C 

C BEGINNING  OF  LOOF  THAT  OOES  IOOITION  ••••* 

C 

JUG-1 

t RCAOIIFILCI  CVIFi 0V2F»0DVFt NLIHO 

BUFFER  INCLFILE.IICOVIF.NLIHOI 
IF  (UNIT  (LF1LC)  «E0»  t)  STOF 

C 

C NEWEN  IS  THE  RECION  FOR  THE  ENISSION  SUNNITION 
C 

S RCADCLFILE I COLDER!  Jl , J-l ,HLIROI 

BUFFFR  IN  flFILEt  II  ICtOEN 1 1)  , OLOEHI NL INO)  I 
IFCimiTCLFILE)  .EO.  II  STOF 
C 

C NEHTR  IS  THE  RECION  FOR  THE  TRMSNXSSION  SUFMATIC* 

C 

BUFFER  IN  CLF1LE* 1 1 I0L0TRC1)  tOLOTRCNLXFO)  I 
IFI UNIT U. FILE!  •EO.  II  STOF 
NWS»1 
6C  TO  361 

SIC  NFNEH  (1I-0L0FNC1 1 * 1 1. O'TRt 1))*NR*0LDTRC 1) 

NENTRI1I»TRC1I*0L0TF«H 
TESTER-TESTER *OVF 
TI-2 
60  JJ-1 

C 
C 

00  111  JFG«1,U 
IF  (JPG  .FO.  Ill  r.O  TO  7 0 
IFCK5S  .EO.  11  60  TO  II 
GO  TO  90 

70  NENEH  < I II  -0L  OEN CRTS)  *C1.0«TR  fill  l*8B*0L0TR  I NISI 
NlNTRmi-TRCXI)FCILOTFCNVSI 
TESTE R-TESTER*OWF 
GO  TO  111 

60  IF  (ENSUE  HI  .CO.  O.lt  GO  TO  81 

CHST0R-CNSAVCC4) 

TRSTOR-TRSAVE  Cl) 

GO  TO  111 

•1  EFSTOR— CLOEN  Cl) 

TNSTOR-OLflTRtl) 

GO  TO  111 

90  ERST0R-0L0ENCNVS-1> 

TRITON  OtOTRCNWS-il 

0 

0 IHTERFOlATF  FOR  THE  010  ENISSION 
C 

110  OLOEHI- 61 C JJ I *EHST0Rf  A2 C JJI -OLDCNC N»S) *63 C JJ) *0L0ERCNV5*1) ♦ 
1 6 1 1 J)  *0L  OEN  C N ¥ S ♦ El 

C 

C XNTERFOlATF  rOR  THE  010  TRANSFUSION 
C 

OLOTAil-Al  tJJI*TRST0R»A2C  JJ)®OL0TRCN03l  *A3C  JJ)*OLOTRCNY5*1) 
1*A|CJ.B*0L0TFCNVS*25 
NtNERCIIl-OLDERl  *C1  ,0*TF  tJX)l*IB*OlOTRI 
REWTRCID"  TRIII1  *OLPTRI 
TESTER-TESTER*  ovf 

111  N»JtMW*l 

IFCMS  .IE.  HUHO-?)  60  TO  110 
EHSAVECll-OlOEfMMVS-11 


012220 

012230 

012210 

012210 

012260 

012270 

012200 

012290 

01230D 

012310 

012320 

012330 

012310 

012310 

012360 

012370 

012309 

012390 

omoc 
012110 
012120 
012130 
012110 
012110 
012160 
012170 
012160 
012199 
012100 
012S10 
012129 
T12S30 
012110 
012510 
012560 
012570 
012580 
C12591 
012600 
012610 
01262. 
0126  3, 
01254: 
012650 
912660 
012670 
012610 
0126  90 

012700 
012710 
012720 
012730 
C127  40 
012750 
01276r 
012770 
012760 
012790 
012600 
012610 
0120  20 
012030 
0126  40 
012610 
012660 
012670 
012660 
012610 
012900 
012910 
012920 
012930 
012910 
012910 
012960 
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>4tnm«io«n amum 


IHIS  PACTS  is  BEST  QUALITY  PRACTICABLE 
COPY  PUR2iISHi£D  IX)  ppQ  . 


ouBWUTtM*  i ho*  tr/tr  opt«i 


ftn  it.Mui  imwn 


: i 


F 


c 


1«« 


u* 


lit 


lti 


IT# 


185 


100 


lit 


190 


19* 


?5B 


*05 


no 


m 


8*1 


twSmt3>»M.DIH<HVS*ll 
EN**ve(ki"ouori(iMvs-ti 
TRSIVEUt-OlOTRlNVS  -It 
TRSIVEftlaOLOTRtNVSI 
TRI»ft<SI»OtDTRlNVSkll 
TRSIVF  < k)»OLDTR  (NVS*8I 

t tcantipuci  ovip.ovep.oovp.nlino 

S IFlFOFlLFILtlt  1*0,180 

BUFFER  lNILFIUC .11  tOVtP.NLIHOI 

irtuHTTurTten  ito.iot.ite 

ttfl  NLIHO"NLtNO*S 
0 Re*t>(LriLEMOLOENIUl,.J»6,NHHOI 

RUFTER  IN  (ItKtfl)  I OLOEN <6 1 , OLOEN{  NL  INftl » 
iMumianutn  i30.ieo.i30 
use  RC*0<LftLeM0L0TR(Jl,J"6,NLI“01 
150  BUFFER  lN<LFILE.lH0LCTR(6t  ,Ot  DTR 1NLIM01 1 

imumxt  tifiutt  .E»).  et  stor 
QucfM<n«cHS*ve<D 
OLCtH<2t»EM3kVttII 
OLOFR(JI"CMSIVf  tSt 

ouctr<h«trs  event 

OUOTR(21"TRSkVF lit 
OLCtR(3l"TRS*VEISt 
NVS"2 

IkO  II" Ilkl 

JJaJJkJ 

irtii  .ox.  huihi  oo  to  tea 
co  to  iro 

c 

C IN  tor  MIS  OCCURFO  ON  THE  FILE  FOR  THE  PREVIOUS  LITER 

c 

18  0 11*11*1 

E**l*  (CM SITE!  3t  kEHSlVEttt  1*8.* 

TRl"(TPS6VE(3tkTRR*VC  <211*0,5 
ICC  NE6EM(Tll"ENlk  (1.0-TRt|Ill*BB*TRl 

MF«tR(II)»TmiIt"TOl 
II-TItl 

IF1II  .UE.  NLIHI  CO  TO  281 
tice  H*nEiMmnwip,v?».,Dv,NUM 

9 NR ITS (H FILE 1 (NF  MEN!  J1  s .)■  1.  NLIM1 

S MRITCtMEUrt  <NEMT«(  Jt  ,J"1,NLI"I 

c 

C NRtTc  OUTPUT  ETLE 

C 

1*8  BUFFFR  OUTINEILC  ill  (V1P.NLINI 
I*  (UNIT  IHEIU  Et  .eo.  et  STOP 
BUFFE*  CUTINF1UE*11  (KEREN  1 11  .NENCN  1NLIN 1 1 
I F (UNIT  IMFIUEt  .EO.  81  STOP 
RUFF*!  OUTINFIUE.lt  (NEMT»(1>,nENTR(NUIN11 
IF  (UNIT  (NFIUEt  .*0.  et  STOP 
IF  ( 6PTS  .UE.  ft  CO  TO  T3J 
00  329  J"1  »NPTS 
VP"V1P*FL0IT(J-11  *OVP 
380  PRINT  908.  J.VP.NFRFNtJt  .NERTRWt 
JENO"NUIN-NPTS*l 
00  338  J«JFND,NLIM 

vp»yip*Fuoi,(  (j-i!*ovp 

338  PRINT  Rod,  J,  VP,  NEWFN(J»  ,NCRTR(UI 
906  PC  UNIT  < 110  »Fl?,  6 ,3C19»  81 
53?  CONTINUE 

6C  TO  <569, 9581  .JOG 
9k*  IRNt"IFNt*l 

IPtlKkL  .LC.  RNUt  CO  TO  160 
SO  TO  tie 
ISO  JU6*2 

508  BUFFER  INtKPlLE.lt (VtP.NUIHt 
IEIUNIT tKFILEI  .CO.  et  STOP 
BUFFER  INtKFILF.  11 1 TR  (11  ,TR(NUIHtt 

iFtuNiTOcrruei  .eo.  ot  stop 

c 

C COMPUTE  THE  BLICF  BOOT  EVERT  1 NIVC  • BEGINNING  IT  VIP 


D12970 
0189*0 
018990 
013CJ0 
013010 
013080 
0130  30 
013:  kO 
01 3090 
013060 
013080 
0130  00 
013090 
013100 
013110 
91 3110 
013130 
eisiuo 
0131*0 
01316; 

cisirc 

0131*0 

01*190 

oiJoor 

J1321C 
013880 
013230 
013260 
013250 
013260 
01 32  TO 
0132*0 
813890 
81  33  DD 
0U310 
C1338C 
01  3330 
0 1 ? Jk.' 
;i3?*r 
C133FC 
0133TO 
C133*r 
813390 
01  3k  00 
JlSklC 
813k20 
013630 
OlJkkM 
013658 
C1S66C 
013670 
0136*0 
Cl  369. 
013*00 
013510 
013520 
013530 
013560 
013550 
013960 
013570 
0135*0 
0139R0 
013600 
013610 
013620 
013630 
013660 
8136*0 
013660 
01 367C 
0136*0 
013600 


i 

\ 

l 


4 


S 

i 
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* 

!■* 


Ibis  pass  is  best  qualiw  mcncABXI 

rouM«tfi?uH»sia»ioaftQ  — - 


SMftOUTlNt  (Ml*  TV  TV  07T*i  TTN  *.6**>0 

I3/1V/T0 

»! 

(70 

stem*  vie 
•OCMIGIN 

r«* 

•oe  us  < 1 1 MtoDv  i ace  xm  i 

i*l*t 

rra  .ct.  it  it  stm> 

01370® 
013710 
013720 
013730 
il 37 VI 
0137*0 

IIUO«IK»»l.l 

0137*0 

m 

xeierciM  .it.  vtMi.oi  go  to  t?t 

01 37 7 J 

1MX  .20.  2)  IOOXCSI 21*  BOOK  Sill 

013700 

CO**OM«l.VV»Tie/T*Vt 

*137*0 

1FICOHON  .12.  *.605)  GO  TO  >10 

013000 

D2NOM«l, 

0130 1C 

GC  TO  2*0 

013020 

toe 

OtHOHal.O'-tXWt-COHONt 

013030 

200 

OVNNF*V 1**1. 02»*/60SI 3. 0” 1 CONON/OENON1 1 

0136V0 

ccr 

201HT  • ( IBOOTrS  It  HO 1 .lU* >1.1011 

0130*0 

#•■00012*111 

013000 

2V  0 

TtSTIBwVlP 

013070 

CH0KG(«*12vO«Wr 

0130*0 

l TEST* 1 

0130*0 

OC  2*0  J»l,MtlH 

013*00 

Tt  J»T*<  J) 

013*10 

IV* 

T*(JI»1. 

113*20 

IF  (TNJ  .IE.  0.  1 GO  TO  260 

013*30 

IF  (TFJ  .IT.  17V. 0 1 GO  TO  24* 

013**0 

TMJ)aCl»HIN 

013**0 

GO  TO  200 

013*00 

2*0 

2*0 

T»«JI«2X2I-T*J) 

013*70 

260 

CONTINUE 

31  3*00 

GO  TO  13*0,350,5311  ,JUG 

013**0 

350 

U»1 

01Vy.CC 

170 

IF  1TE5TE*  .IT,  CH1NGEI  GO  TO  1*0 

CiVPIO 

21* 

COHON»CH»NGE-BOP 

01*070 

B1»#OOIES  U.TEST  1 - (000  IF*  UTEST  > -BOOIC  S UTE  ST  *1  > l*COHON 

01*0  30 

CH*NG»«CMANGE*DV*er 

ji*jv; 

IF  1 CHANGE  ,CE.  BOOM.0I  GO  TO  30C 

C1*0*0 

GO  TO  1*0 

oi*:  6; 

2*0 

300 

LT*ST»ITF*T*1 

oi*. 7: 

«O0»0OOM.  0 

SlV.05 

1*0 

CONTINUE 

01*0*0 

NVS»N*1-1 

01*100 

c 

•1*110 

2** 

GO  TO  60 

01*120 

210 

OEM  IMP  XFT1E 

oivi3; 

c 

C6U.  StCOMPITIMFl) 

01*1*0 

;i*l*0 

TI«»TI2tl-TIVE 

01*160 

2»0 

»RIMT  7*0,TTN»1,TIH 

01*170 

HETIMN 

01*100 

c 

01*1*0 

c 

01*200 

c 

01*210 

27* 

220 

FO*N»T<*0  THE  TINE  »T  THE  SHOT  OF  ENUP  IS  *,F12.3> 

01*226 

2ve 

F02NATI*®  TMt  TIME  »T  THE  END  OF  ENUP  IS  *,711.3/212.3. 

01*230 

i*  secs  vtur  vtoumo  fo*  this  addition*) 

01*2*0 

c 

(NO 

31*140 

01*260 
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-u.^. 


W]tn " ViMmi1 1 1 m 1 1 


IBIS  MU®  IS  MP*  'I'SilCW© &tM 

OOPX  fURWSHSD  10  >i*Q  — — « 


tttotvttnt  fwttw  ?*/r*  o*t-i 


tr*  list*/?* 


t 


i 


it 


it 


10 


t» 


jo 


w 


00 


Of 


to 


» 


00 


** 


7# 


r* 


SUttOUTlNf  IN90ON<ntTS'|TVtC.lfUtt*mt> 

COttOn  CHOP*  (tot 01  «010T*  ( 1610  ) «Mt  Mint  0*10) 

cotton/nop*/  umi.iotMiiOittiook' 

conNon/xtnm/  NNa.,xiot7>,jrciNT,*m,Y*vr,»*otjo(T>,*((T», 

C DV)«l,«l,OUUO 

cotton  /o nm#  ovit,OYtt*ooY*,niio 

cotton  /ntMt/  *it,**t,o»t,Niin 

DttCntlOM  t*J4Ytl*>.*i<l#».*»<lHt»3tiO>»»*<l(?,TM»VII*> 
01NWSI0N  T»(»*l*>,)«WT*(l>,*OOItJ<l»l> 
tOl'lOOltMCt  <TftUi»N«MTMl>l 
MU  MrNCNt<Knr« 

c 

c 

o •*•••▼►»  iut»ouTiMt  cess  tns  ut«  oooinon  too  tnt  mission  »**••*• 
0 ••"H.OOOttO  toon  snoot  to  MOUM  •••••*• 

r. 

Cxtr<Tl*fXt.*M4xi(»ir*.i,Y>> 

»ionT(Mi*(i.i»tet-it****J)/(T»tra.i3TM»x/T»vt)'i.O) 

lx*Nin»tx*(-i?*.  t» 

cou  sceo*om*n 

mint  tio.tint 

l*ni»t 

mint  this  output  is  numen  on  out  “.tout 

c 

mint  *,  - moo***,  true.  iftic.Krtit  “♦•‘flit  ,i  tut  ,«mr 
nni*«*o»ti 
10  otnino  nrm 
nmiNC  iriu 
ncmnc  true 

•urn*  iNii»ne,t>  ixiotn .*xom> 
it  (Ukimruti  .tn.  t » sTot 
nuttet  iNiKtiiFtiMxxoat.nioveni 
ltiunitmtiict  .e*.  0)  Slot 

ewuntiit.ii  itnut  »Ni«vtt> 
xtiuwn  (nttiti  .to.  oi  srot 
00  JO  K«i,o 
T *3  Mt  (*>•(. 

so  rn,uvtt«»»o. 

4TYtl*I Tttt 

TMITYtt  01  60  TO  030 

C 

C 1/1  MTIO  ONLY 
C 

JU6»S 
J06»f 
60  TO  3«0 

nutrto  iNante.n  <cvtt,nu*c> 
rttiniTarrici  .to.  o>  srot 
■ ctrr*  JMdtriE,  lMCtO'HC  t»  .oiocntnirtou 

IMUNIT  (17110  .to.  0)  STOt 

store#  Tnatuctt)  (oicTnit)  fOiomiNiito) > 

X#  (UNIT  (It  II  Cl  .C  ..  0»  STOt 

to  ojo  i«i,ntxn 

£nrnu>«T«txi*oioct«i»u.B-Tn(i»»*eo 

:nTtrx)«Y*(ii*oioT»n> 

TESTft*  TCOTCtOV* 

Tt (TtJTtt  .IT.  CHINCCI  60  TO  930 
ON'ONaYHOnSt-YOD 

NtnOOCTCStlTCSTi *1000 IfS  llTKTI-BOOieS  (lT£$T*l ) »*CO*CN 

eHtnc(«CH«not*ovoor 

xntcniner  .tr.  too«t.ot  go  to  too 

to  to  sso 

*40  ITf JT.ITIST*! 

BQ0*  000*1.0 
oso  :>  nTim>t 

PC  TO  100 

oft  x«ni»wni*i 

IMTMl  .It.  XNU  60  TO  JOB 
GO  TO  tie 
c 

C 0U  MTIOS  rxet*?  1/1 

c 

030  ll*ITVM*l 
J06«l 


mm 

oiotos 

Hitte 

01OJ01 

*10310 

(lost* 

*10)1* 

*1030* 

010310 

61030* 

(lost* 

*103*1 

6103tt 

*10**0 

*1001* 

tiooto 

010*3# 

*1000* 

*10003 

*1000* 

*1*0/0 

0100*0 

#1*00* 

Cl**## 

C1O010 

tiooto 

810*30 
#10*00 
#10**0 
#10**0 
not  rt 

#10*40 

#10*1# 

*1**03 

tl*61t 

oi**re 

01*630 

*1***0 

01**90 

01**** 

;i**t: 

tl*b«0 

01*6*0 

01*70. 

01*710 

01*7tC 

01*730 

01*7*0 

01*790 

01*760 

01*770 

>,1*7*0 

0107*0 

Cl*#>,C 

C1>#1* 

010770 

ei*M0 

cio#*. 

01*#*0 

01*800 

01*070 

01*480 

01*4(0 

01**00 

(1**10 

oi**ro 

01**30 
0 1***0 
01*0*0 
01*9oG 
01*970 
01**46 
(1***0 
019000 
019610 
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ans  PACK  18  9KST  QUA UtTt  FSfcCTiCA*» 

UWMOOfYfURMSHSDlCDOC  


MWUTtKI  (MBC«M  tvn 


»tn  *.**tt i ii/i*/7i 


•• 

•« 

*• 

111 

lit 

111 

lit 

uo 

ltt 

191 

13» 

It* 

ltt 

It* 


ttai.i/imn«i.ii 

00  *•  J**M,XTtlt 

IM*JH 

Ml.HW'tHI 

POINT  •«”  ►•  “,P 

e iTirt  is  T»»f  nunioiton  or  thc  mtxo  or  oiooy  to  Ntto*  nmicn  is  closist 

C TO  * 1/1,1/l.S.l  (TO  TO  Cl/1,  r tOUlUS  TNC  NTClMoClL  01  TNXS  NITXC. 

C 

c ••••»••  thi  rot  ion  in*  m tmc  const  ints  ro*  tk  imunii  * roiNT  t 
C IITFWOUTIM 

UIJH)  MMIH.Il'IM.II/l.l 

«ww»(r**H.iiMH.ii*e.i 
*sijM»»«r*tru.  •»•»»-*.  •»•*.« 

it 'jtt>«t*(»**(-i. •>/*.* 
to  CCNTIWt 
c 

C •**»•*•*•••*  IIIIWIIIC  or  too*  THIT  Does  lOOITXOk  ••**•*•••••*••••* 

c 

J06«l 

S ftticttrut)  0»ir,0»fr,00»r,*«LTN0 

•urrr*  XNiiriuMHCVir,  nlxmoi 
irumiT  anti*  .to.  it  stqo 
c 

C TtlNSMSSION  IT  Tttt  MtStNT  UTtt 
C 

s MticarxuHOt.ocMM>,j«i,«iLrio> 

tn  turn*  XNarxu.iHOtocNtn.oiorNtNLXNou 
irtuMTTarxtti  .to.  it  sror 
ourm  iNiirnt.1  MOtoTi in  .otoTt  iNtinoii 
irtUNirarxtti  .to.  «)  STor 

sti  NVS'l 

60  TC  ?*• 

3*0  NE«CN<l)aTIMi;'»Ot(iCI><l)M1.0*TI<l»>*Bt 
NMTMS  »«TOI  1)  *OLOTI  1 1) 

TtSTtttoTtSTftMlwr 

xrn 

M JJM 

C 

c 

00  1*0  vP6M.ll 

xrijrc  .to.  U>  60  to  TO 

imnvs  .to.  it  to  to  *o 

60  TO  •* 

7T  NtNtNIXI>»TI»m>*OtOtN<Nn>t«1.0-T*<Il>>»** 

NCNTNII1I»TNIXI>  •OtCTUNISl 

TfSTrO«TtSTt»«Otr 

60  TO  111 

•o  xrirMsittiti  .to.  i.o»  to  to  «i 

f«STON»tNS**tl*l 

TMTOt»T»SI»t«t» 

60  TO  110 
II  CNSTON»OtOCNIt» 

TMT0H»0tOT^HI 

to  to  no 

tr  trstoi»c«.otmi«vs*i> 

TN3TOt*OtOTNINVS“  j,l 

0 

C INTttrOllTICN  Or  THt  NtlINT  TNINSNXSSXON 

111  OUWIUl(JJIMMTO**«IJJ*,OlOlHO«XH»l(JJ>,OtOIHtW5»lM 

J 0l0TtT»61  l5r*TISTOI*IK  J4>  *OCOTMNVS>  tlSl  «W>  •OtOT*IN¥S*l»  ♦ 
till  Ml  •OtOTt  INIS*  t) 

•»tttiMXx>»ti.O“Tttxxn»tt*T«m»*xot»ii 
NtMTIMll)  •VtlXX»*XOT»X 
«STtS»TtSTt*tO¥t 
Ut  HTtMrttM 

xrm*t  At.  NtXWO-O  tO  TO  ltl 
(NSIVttl  V-OtKN<  HTt*t> 

(Msivtm-xotHiNvsi 

rNS»vcts>»ototHi*m*i> 

tNSI»ttt>*ttO(N(NVS*tt 

TM*VHt>«CtOT*INYS-l» 

TMIttltlaXOTIMNVSI 

TtSI*«IJ»«XOTMMVS*U 


1110(0 
titoxo 
ouite 
imit 
Oil CIO 

oitore 

91 toio 

me  to 
01*110 
01*110 
01*1(0 
01*110 
01*110 
01*1*0 
01*110 
uiiro 
01*11! 
oitit; 
01  (too 
omio 
ciscto 
01*00 
01*110 
UKIC 
01*00 
attire 

01*09 

01*00 

01*133 

cum 

01*1X0 

oitu: 
01*1*0 
01*1*0 
ei*oo 
01*17. 
01*1 19 

ci*i*e 

3t*fcvV 

ji»*i: 

:i**x. 

Cl**lt 

9i»**r 

vl**S. 

oiotir 

01**70 

11*110 

llttto 

01*900 

31**1! 

n*»!i 

01**10 

•1***0 

11**93 

unu 

01**70 
01**00 
?1***C 
01*110 
01**10 
01*1X0 
01**10 
01*1*0 
01*1*0 
31*1*3 
01**70 
11*111 
oiuto 
*1*711 
01»T10 
91*7(3 
01*710 
• 1*7*4 
ll«7t* 
9i»7*J 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

J350M  COPY  FUKNlSHJfiD  TO  DDC 


SUB  SWT  INC  EMMNN  76/26  OPT«l 


ptn  stttti  ii/iwri 


ill 


i*o 


i*s 


in 


1M 


III 


US 


1*1 


111 


no 


m 


m 


in 


no 


TftSOV(t6)aULOT*<NVS-t) 
t ttoPuriiei  wiP.wor.oove.NLtNo 

t iMtwunmi  iMi if* 

iwii  xnilfilc,imovip»nlxmo> 

SMimnartini  in, in, no 

in  nlipo^nlihoo 

I M»CUFHIH0L0fPIJI,.P'6,NLlH01 

•UPPtN  XNa»U.C,l>(OL3IM(6l,OtOCH<HLXNOIt 

XPiutmariLtii  imimii 

II  WMII  XNUFXLt.tUOiBTNU)  .OLOTMNLINOII 
IKMITMntll  IMIMII 

63  0L01Nll)aCHS«VC<ll 

OLOtNItlaCNSOVCItl 
0L.0EMt3)aCNM««SI 
OLCTNtllaTBSlVf  111 

euimti*  mmtn 

OLOTNIStaTMOVCt*, 

on*i 

1,60  XX>11«1 

xrixx  »6T • NLIM  60  TO  too 
60  TO  ire 
c 

c on  cor  has  occuteo  on  the  rxu  roe  the  envious  loye» 

c 

100  II«II*1 

eni*<chs9vc  nmnsovE  itn*o.9 

TP1«  (TPS* VI  <3>.T6S»Vttf  >1*0.9 

no  NC«N«u>»u.o-T*iixn»et*TNUi>*ni 

NENTPIH».TN1*T9IIII 

TlaXI*! 

IPlXl  .It.  NLIN1  60  TO  (00 

e 

C NNXTE  OUTPUT  PILE 
C 

mo  NPITEf HFXLC)  V1P,V2P,CV,NL1" 
t M«ITMHmn(NEMfHMI,>l,NLXN) 

t mmiNPuri  iNfNTNui  ,j»i,ni.ih) 

190  BUFPf * CUTIHriLF.lHVlP.NLIPI 
IFIUNIHPFUE1  .10.  01  STOP 
BUFPEN  QUT'NPILF,!)  <NEMEH<l).NEP'*<Nt.XH)> 

IPIUNTT  (HFXLC)  .CO.  01  STOP 

9UFFE»  OUT  (NFXLC  ,1 1 INCPTM1)  ,NCPT*(NLXH>> 

XPIUNIT  WPILC)  .10.  01  STOP 
IPXVtaS 

XPINLXH  .CO.  91  XFX«P*1 
IP  I NPTS  ,LC.  0 I 60  TO  333 
00  ns  Jal.NPTS 
VP.V1PPPL0»T|J-1I*CVP 
Sto  POINT  900,  J,VP,NL'NCN|JI  tNCNTPUl 
JFNP*Nt  TN-HPTS*  1 
IPIJCKD  .IT.  0)  JEN  C*NLT  9 
00  3X0  J»JENO,NlIN 
V*»  VIP*  MOOT  0-11  *DVP 
ISO  PNINT  900,  J »VP ,NW IN  (Jl  ,NEPTN  (J ) 

900  PO9P«TIX10,Plt.*,3PlS.SI 
333  CONTINUE 

60  TO  1969,990)  ,<106 
969  XKM.9lPNl.Pl 

IFtIKMt  .LE.  XNLI  60  TO  190 
60  TO  CIO 
190  JU9»2 

390  BUPPCP  INIKPILE.l)  C V1»,NLIN) 
tPtUNIT  IPPILC)  .EO.  0)  STOP 
BUFFI*  INIPPItE,!)  ITNIl)  .TNINIIN)) 

XPIUNXTIKPILEI  .00.  01  STOP 

C 

C CONPUTE  THE  BLOCK  BOOT  CVCPf  1 POPE  • 9E6INNIN6  OT  VI* 
C 

BE6IN9V1P 

BCB99ESIN 

I9l 

ore  pooiesi T) »bboot  ibecxni 
s»t  *t 


otirro 

OlSTK 
019190 
919009 
019016 
0190 tl 
019030 
119060 
019090 
I1S990 
1190. » 

noon 

(19090 

J10900 

019916 

019920 

119910 

119960 

*19990 

.19960 

min 

3199IC 
019990 
019000 
31901C 
3196  20 
019030 
319360 
016090 
019060 
C160T0 
31930C 
016090 
316100 
016110 
.16120 
016130 
C1616T 
316193 
016160 
:i6ir: 
01 61  60 
016190 
016200 
316210 
016220 
31623C 
11626. 
.16290 
r 6260 

u 1 62  TC 
016283 
016290 
016300 
016310 
016320 
016330 
016363 
016390 
016360 
0163T0 
016300 
016390 
316630 
016610 
016620 
016630 
016660 
016690 
016660 
0166T0 
016600 
01*690 
.19930 


•mm* ipw***'  til  r -■*» ■-  >4 '•'WSVt  i"?,|!'  •M«r>i:S'-,-Wi''1 


«^*r‘****'V'*'v 


tills  PAGE  is  B®S«?  QUALlrrmCTICiBMI 
ISOM  COPX  PUKttlSHM)  XODD.Q  — 


SIMMWtlMI  I NDONN 


r»/t»  o»T»i 


ftn  4.**4to  itm/n 


irr 


too 


tso 


IM 


MO 


f33 


1*0 


too 


tn 


tto 


rrn  ,6t.  m>  «t» 

innt>»niN«i.i 

IMMMN  .IE.  OtMl.ll  «0  TO  >71 

i7oi  ,ro.  to  •ootctot»*eooi(iut 
ccc  ppioit  *»  nooxit  ox 

eOHCN»1.44‘M17/T*Vt 

IFICONON  .11.  4.009)  00  TO  >00 

DtNON>1.0 

•o  to  toe 

toe  DtNON»i.o*E»77i-coHON) 

too  0 *007«V»M1.  0E-4/M3 0 S. 0 - i COHON/OINON)  > 

OO*  0001  tStl) 

TXJTH»«* 

CHOMO  (*01000007 
LTttTal 

OC  MO  >l,NttH 
T7J»T«(J) 

TN?  J>  *1. 

X7  « TOJ  .It.  0.  I CO  TO  t*C 
X7  ? T7J  .cT.  174.0  > 80  Tn  t*0 
TOtJtaCXOHXN 
00  TO  ttO 

tOO  TN?J»«E*7(-TN.» 
t«0  CONTINUE 

37C  rCNMTI  1X1710X11. 7) 

00  TO  (740.7O9.fS8l  .JUG 
300  11*1 

170  I7(TESTfO  .IT.  CHONG*)  GO  TO  HO 
CCH0N»CN*N6E-70P 

BMtOnifS  (ITFSTI  - (HOOT? S (UTEST)  • BOniES  (LTEST  ♦!) ) •OOMON 
r NO  NGf*CH*NGf  ♦UV487 
IF  (OUNCE  .Gt.  00041. C>  GO  TO  300 
GO  TO  100 
300  (.TEST*  TEST*  1 
OOOrROOU.O 
100  CONTINUE 
NVS*N0?-1 
GO  to  00 
210  worn  H7TLE 
c 

COLL  SrCONO(TIHfl) 

TXN-TINfl-TIHE 
77INT  7*0  » TIMtl .T JH 
METU7N 
C 
C 

ttO  70NNOT ( *9  THE  TTHE  »T  THE  STiNT  OF  tNDONN  XC  *.Fll.3> 

t*P  70NH*T(*0  THE  TIHF  »T  THE  ENO  07  EHOOMN  XS  • » Fit.  3/Flt .3  . 

1*  SECS  HETf  NEOUXNEO  FOR  THIS  400ITI0N*) 

5 

ENC 


01*010 

0101  to 

01*030 

01*040 

01*080 

ti*o*e 

01  >070 
01*000 
oitooo 
01**00 
01**10 
omto 

01**30 
010*40 
L16630 
01**40 
01»*»0 
01**00 
oi**oc 
01*730 
31*713 
Pit?  to 
01*730 
01*743 
01*730 
01*760 
C1*??C 

c;*?30 
otfcTor 
GltbOO 
;i*ci‘ 
C16»?l 
01*730 
01**60 
01*030 
01*003 
01**7  3 

oiteot 

tl***' 
01*900 
316910 
3169 tO 
016930 
01*940 
016030 
01*9*0 
01*970 
01*9*0 
81*943 
017.03 
017010 
0170t0 
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137 


5F 


V 


IBIS  PA01  IS  BSST  QUALITY  mmCABU 
ggpjt  COPY  PURJilSHKU  TO  DDC 


SUBROUTINE  MlRACO  T4/T4  OPTpI 


pth  ».»«»(•  n/u/ri 


f 


i 


s 


it 


70 


H 


30 


31 


10 


49 


50 


55 


50 


55 


TO 


SUBROUTINE  MlRACO 
C 
C 


CAL  COLITIS  MONOCHROMATIC  ABSORPTION  CQEFFICmT  FOR  SINGLE  LAYER 


• USES  HXRACC  ALCORXTHH  FOR  THE  OORPLCR  PROFILE  • 

* • 

C 

e 

COMMON  CMU<*5S) , S <790 ),  ALFAO  <7503 ,EPP<750>  ,P0l<75B>  , 

C CFOPTHIf5»».»ErAL»IE50» 

COMMON  FF |S500> 

COMMON/  MAIN/  FFILEi  KPAHELt  P0  * TEMPO 

COMMON/ NfN/  NMOL*  XIO(T)  ,S£CANT,PAVE,TA Vf  • MMOLID •. 7) ,MK < 7)  , 

C OV* Vi* Vf*NLATER 

COMMON  /XSUB/LXMIH,XLO*XHX  ,VtOT*VTOP ,VFT,I EOF ,IP ANEL  ,XSTOP,IOAT A 
COMMON  /NXX/  NNF,OXF,NF,NWS,DXS,NS,NMVS,OXVS,NVS,NFMAXtNSMAX,NVSH 
1AX 

COMMON  /SUM!/  MA  XF*HA  XS*MAXVS  *NL  IMF  *NL  IMS*  ML  IN  VS  * NLO*  MHX*  OVS  * OVVS 

COMMON  /XPANCL/  VlF, VIP*  0VP,NLXM, NSMIFT ,NPTS 

COMMON/XT I ME /T INI *T  XMROF , TXMCNV • T XMPNL 

DIMENSION  ALFCORIT)  .ALFOltri  *U(7)  ,SC0R<7) 

qiMENSICN  *0<751> 

C 

DATA  NNF  / 4/*  OXF  / g,f;/,  NF  /t01/r  NFMAX  /E51/ 

DATA  NHS  /15 / , OXS  / 0.00/.  NS  /EOt/i  NSNAX  /Mi/ 

OATA  NHVS/54 / . OXVS/  0.??/.  NVS/701/,  NVSMAX/M1/ 

C 

DATA  I ENTER/0/,  LINIK/750/,  NSHXFT/S7/,  NLINF/2401/  * NFTS/  0/ 

C 

OATA  SUB1D  /IBM  HIPACO  / 

C 

IF  <IfNTEr.NE.0l  GO  TO  10 
TENTER*  1 

NLIMS*<NLI“F/4>*1 
NLIMVS*  INLIMS/4)*! 

C NOTE  Irxvs/OXF)  IS  15  AND  <OXVS/OXS>  IS  4 
NBOLNO*  EOXVS/OXF)  •FLOAT  INMVSI 
MAXFaKLIMF.MMCUNO 
, MAX«»<MAXF/NI*i 
HAXVS«<MAXS/4)«1 
CALL  SHAPEOIXOI 

CALL  MOLECCt.MMOLIO, NMOL, TEMPO, TAVE.PO, FIVE. SCOR.ALFCOR.ALFOi) 

10  PRINT  500 

PRINT  R?0,  SUB  ID 
T IMPDF*TIMCNV«TIHPNl«0. 

PRINT  544,  (XIOII)  ,T»1,T> 

RE  NINE  KFU.E 
IECF-0 

PRINT  510,  KFILE 
PRINT  515,  SECANT 
PRINT  515,  PAVE, TAVE 

PRINT  53f,  <HMOLID<M<  ,NMN>  ,M»l,NMOl» 

PRINT  540,  *V,V1,V7 
OV»»OV 

OVS«<CXS/OXF> *DV 
OVVS»<OXVS/OXFt»OV 
BOUNOmFLOAT  < NBOUNOI *0V/7 . 

80 U HOF*  BOUND*  <nXF/OVVSI  /?. 

PRINT  »4f,  ROUNOF 
ALFHAX*  IJOHO/FLOAT  (NNF) 

Nl 0*NB0  UNO*l 
NM I«NLI HF ♦NSMIFT- 1 
DO  90  1*1,  MA  XF 
50  FF<I)«f. 

C 

I WRITE  < KFILE)  < XIO<  I)  ,1*  1, »)  .SECANT , PAVE, T AVE  , <HMOL IO<  W »N»t  ,7)  , 

* T <W«<I)tI-l,NMOL>,DV,Vl,VE 

BUFFER  CUT«FILr,i)  <XID<1»  ,NLAYER» 


OIBIOB 

00(110 

0001(0 

100131 

050140 

001150 

000150 

000170 

oooiai 

000150 

tcotai 
000710 
0007*1 
000730 
000740 
000750 
0CC75B 
0Cu770 
0007B0 
6C  0750 
000300 
000310 
0C0370 
3CC330 
000340 
0003S0 
Of  0350 
0(0370 
000300 
0(0350 
0(0400 
000410 
000<70 
0004  30 
C00440 
000450 
00046C 
0C  01,70 
0(0450 
000490 

ocosj: 

000510 

000570 

0C0SJC 

0(0540 
000550 
000550 
000570 
1(0550 
000540 
00(500 
0(0513 
000570 
OC  0630 
0CC640 
OC  0650 
000660 
000670 
0306  OC 
00069C 

or  0700 

000710 
001770 
005730 
000760 
000750 
005760 
0(0770 
100700 
0(0700 
MMOl 
001010 
It  OAK 
1(0030 


n 
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JB/M  COPI  FUKWSHBD  IQICQ  •— p-" 


tuavourttt  Hiatce  tn/fn  oft.i 


ftn  n.unei  iim/n 


TS  I*  (UNJTfKFXLEl  .CO.O)  ST  OF  “W*  (MM 

C KMU 

VFT»A 1“ f « • BOUND  0008(0 

VB0T»¥1 -POUND  DC  08  70 

VTOF»m*»OUND  DC  0(8  0 

•0  C 00  08  90 

XKTOaO.  F99l*TCMP|  000900 

xkt  •(.»m,T«vr  5ic»ic 

XKTFAC>(1./XKT>  - tt./XXTOJ  000920 

e«U  90UC<2iHN0U0,NP0L.TENHe,TAYE,  FC, FAYE,  SCO*, ALFCOA,  At  FDD  000930 

•9  00  M NpI.NNOL  0(0990 

00  U<N»«N»UN»*$eo*  IN) ‘SECANT  0(0990 

ICNT-0  000990 

SUNAtFM.  000970 

NCHN6«t  00  0900 

90  C 000990 

90  CONTINUE  00 10 00 

c come 

100  CONTINUE  01 1C  EC 

C 001530 

99  CALI  SECOND  (TINfO)  001040 

IF  <IEOF.NE.0t  60  TO  ISO  001090 

CALL  AOFILEIGNU.S.AlFAO.SFF.NOO  0010(0 

CALL  SFCONP  <TI«f»  001CT0 

TINPOF*  TIHAOF.TIPE-TINEO  0C10SC 

100  C 001090 

TFdfflF.KC.OI  60  TO  1(0  0C1100 

c nemo 

C NOOIFT  UNE  DATA  FOF  TEMPEAATUAE , PAESSUAE , 8NC  COLUPA  DENSITY  OtllEC 

C 001130 

109  00  170  I«lLO.XMI  001140 

H»90L<tl  001150 

»FOPTmn*su»»u«*o  qciko 

IF  tEFDPTNd  I ,LE«  0 • I 60  TO  139  0P117C 

ICNT»  ICNT  *1  00 1 1 80 

110  ALFO»GNUII)*AlF01(H|  001190 

TFTAtF0.6S.09l  GO  TO  110  5C123C 

»9I»'T  91.9,  GNUID.S  Ttl.Uf  APU)  .AtFD.PV.H  001E1P 

•LFO«OV  0 C 12  20 

NCHW.*NCHNGM  OC  1230 

115  110  IF  (ALFQ.tF.ALFN(x)  CC  TO  1E0  0C124O 

P*INT  550,  6NU(I1  ,STI),ALFA9tI)  ,AtFO,ALFH»X,<*  C012S0 

AVFn.AIFNAX  OC 1260 

NCNNG*NCMNG»1  0C12T. 

1E0  CONTIUW  CC12  30 

120  5UHALF»SUPAC  FyALFO  001290 

NSCALFUl  •1,/AtFP  0(1300 

EF0PTN(I)»EF0PTN»I>*EXP(-ffF(I>»X9TFAC>  ‘AECAlF(I)  • 001310 

* U.-CXP»-GNUtI)/TKT>>/'<l.-tXP(-GNUtI>/XKTO>>  OC13EO 

ISO  CONTINUE  0(133; 

129  IF  (NCHNG.GT.1001  GC  TO  19"  0013*0 

INC  CONTINUE  OC  1350 

C 0013(u 

CALL  eONVFNOIGNU.EFOPTH,  XSCUF.POL.FF.XD)  001370 

C 01138  0 

130  IF  ( IPANEt  • EG .01  Go  TO  109  001390 

C 001N90 

190  CALt  P'NEIO  (FF.XFUE)  utINIC 

f QC 1420 

CAtt  SECOND  (TINE)  001430 

139  XPANEl»XPANSL ♦!  00144C 

IF  (ISTCP.NS.l)  GO  TO  1 NO  001450 

END  »ILE  XFItS  001460 

1(0  CONTINUE  001470 

"NIPT  955 ,TXNE .TIHACF ,TINCNV ,TINPNL  0(1480 

140  IF  (ICNT.LE.OI  GO  T C 170  C01490 

A YAtF«SUNAtF /FLOAT ( ICNT)  001500 

FAINT  9(0. AAAtF, ICNT, NCHNG  001510 

0 001520 

170  AETWN  001530 

145  r 001540 

c 001550 

C 0015(0 

990  FOAHAT  UM0»  001570 

904  FO»HATIiX,74tO»  0C1580 

199  910  ro#HATI»P  OUTPUT  FILE  NO.  *»  I5>  01 1590 
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SHIS  PASS  IS  BEST  QUALITY  PRACTICABLE 
FfKU  COPY  FURNISHED  TO  DDC  


JU«  HOOT  IH*  HIPACD  7*/7*  OP  /■! 


PTH  *.t ♦ ***  0l/l*/7* 


19k 


UO 


1*9 


919  P0*H*T<*C  SECANT  *•  719.91  OC16C0 

• 70  POPHAT  UO  (1PE10.3)  I 001610 

979  POPHAT<*0  PPCSSftlS)  PIT. 97  *0  TEhPIK)  »•  711.*)  0016*0 

970  73ANATT7SH0  COLUNN  DENSITY »NOLECULES /CN*>**)  , 001610 

1 / »9X,»6,  • • • lPtia.S)  ) 001660 

960  7QPNATX*.'  OV<CP-i)  ■ *71*.6/*9  Vl(CH-l)  ■ »Pl*.*/»0  V*«CH“1>  ■ * 001690 

1 lit.*)  001**0 

9*t  POP HAT) 17H0  tOUNOP  TCK-l)  * »P0. *>  001*70 

969  POAPAT  C*  — ..................... —•P10. *,  El*.  3,  3710.6, 19)  001600 

990  POPHAT  0*  ♦♦♦♦♦♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•♦♦♦*Pll.*,Ei*.3»3P19.6,X9)  0il*90 

999  POPHAT  UNO, 10X, ‘TIME*, 11X,»PEAD»,*X,*CONVOLUTXON*,IQX,*PANEL*/  001710 

1 *719.3  I 001710 

960  POPHAT  0*0  AVCPAG*  NXOTH  ■ •710.1,  0017*0 

C • NO.  LINKS  ■ *110,  • NO.  CHANGES  ■ *110  ) 00173C 

970  POPHATO 1H0,PX,A10)  0017*0 

CNO  011795 


SU9P0UTIHE  COHV7HO 


7 */7*  OPT«l 


TtN  *.  6*  *76 


11/16/76 


10 


19 


70 


79 


30 


39 


*0 


*9 


70 


TO 

*0 


90 

60 


SU9P0UTINE  CONVPNOlCNU, EPOPTH*  PECAL7*MOL»P7 , X0> 

COHHOH/  HE*/  NHOL  ,*1017),  SECANT ,PAWf  ,TA»E  ,MNOin  C7)  ,MK  0 7) , 

C TV ,91 . V* ,HL  ATEP 

COHHOH  /XSU9/LIHTN,  ILO,  INI * VPOT  ,VTOt  ,V  PT  , IE  07  ,IPA)«  L , 1ST  OP  ,t  CA  TA 
C01HON  /H**/  NM7,OX6,Nr,MMS,0*S,NS,NNVS,DXVS,NVS,N7MAX,NSHAX,NVSN 
1 AX 

COHHOH  /SUBl/  HAXP, HAXS,  HAX VS,NL  IHP ,NLI7., NL I7V9  , NLC  ,N7I  ,OVS,  CWS 
CONMOH/XTIH7/TIHE .T1HP07, TIHCH9,T IHPHt 
OIHfcHSICN  6NUUI  , EPOPTMU),  PECALPU,,  77U) 

OIHrNSIOH  t/OL  1 11 
OT-EHSICH  XOU) 

CALL  S7C0NP  (TTHEO) 

•AT»x«rvvs/oxws 

7VS«7lOAT<NMVS)/1X>'« 

C0N7-CWWS/nv-t. 

17  (ILO.GT.THI)  GO  TO  90 
00  69  I«XlO, INI 
Of  PTHI»E7P*TH(I ) 

I7fOrPTMT.LE.n.)  *0  TO  *9 
7SLOPC«P£C»LPm»PATV» 
ttHTXGNUttl  "W7T)/OVVS 
87VS-WS/7SLUPE 
JHAX*  17IHT  *97  VS ) ♦*, 

IP  (JHAX.LE.HAXVS)  7.0  TO  70 

ILAST«I-l 

60  TO  90 

JHIN.«7IHT-arVS)*t. 

J7IN7»|lIHT»COH7) 

7 VS*  < PL  CAT ( JHIH-?) "IIMT) *7 SLOPE 

IP* TVS* t7LOAT<JHINP)-XXNT»CON7)*ZSLOPE 

00  70  JJ*JHIN,JHAX 

J7«JHIH7*JJ 

7P«77»7  SLOPE 

I7P«A9S  IIP)!  1.  9 

PPTJP )»PP IJP )♦ DEPTH T*  XO  IIEP) 

CONTINUE 

COHTIHUE 

1IOATMOIPOP  NOPE  DATA  PDUIPEO  ••  IOATA»I  IP  HO  HOPE  OAT  A PEOOIPEO 
X*  A 7f  l»  TO  AT  A 
60  TO  60 

continue 
IPANEl«l 
ILO*IL»STM 
CALL  fPCOHO  1TTHE) 

TlHCHV»TIHCHV»TIHE-TI*>FO 

PPTUPN 

PNt' 


DC  1760 
01177C 
001760 
0017*0 
OCIOOO 
OtlSlC 
00167P 
011630 
3111)1 
001690 
OC 16  60 
0L16TC 

ocieer 

ooi»vc 

0C1°  3t 

uc  s 9i r 

001P7C 

PC  19  3. 

Oi.lP*l 

001990 

001960 

0P197C 

601960 

001990 

or  7i> U U 

CP  7010 

007070 

0070  3: 

OC  70*0 

017090 

00706" 

01  7070 

Of  *0  60 

007090 

0071C0 

007110 

00*170 

00713C 

com. 

00*190 
00*160 
00*170 
00*160 
007190 
i(t*oe 
017*10 
Off*  to 
10**30 
Ol*t*C 
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■ar  S3 


■MIS  PAG!  IS  BEST  QtJAI/tT*  PBACTICA&S 
JBO*  COP*  SUKJtaSHKD  IQ  DM 


sumoutinc  MWi.0  79m  opt«i 


PTN  9.6*926  03/19/76 


1 


9 


10 


19 


20 


29 


30 


33 


49 


90 


fuowutinc  pmoclo  cff.kfilci  ocimo 

COMHQN/NCM/  NHOl,Xin(7),3CC6NT,P*VC,TIVC,PN0L10(T),NK(T>,  OC2260 

c ov.n,  w,w.*Tt»  002270 

fONMON  /VSUI/LIHIN,2LC,IHI,V90T,VTO*,VF7.ieOF,IPftNCL,ISTOP,ICtT6  002260 

coupon  /mi/  H*xF,wivs,N*vvs,NtiNF,NLiHi,NLiFvs»Nto,NHi,ovs,ovvs  002290 
COHP0N  /IPANFL/  V1P,V2P,0VP,NLIH,NSMIFT,NPTS  002309 

CONNON/  VTIHC/TIMf  ,T  XNOOF  ,TINGNV»T  XHPHL  002310 

DINCN3ION  FP(1»  Ok 2320 

cm  sccono  (tincoi  002330 

ISTOP«0  002390 

NNHIa(V2-VTT»/DV*1.9  002390 

IF  (NHI.6C.H9Hn  IST0M1  0C2360 

IF  (ISTOP.tQ.ll  NHXaNWI  002370 

JHXP»«NLO*NPTS  002300 

JJPP«NNI-NPTS  0C2390 

LlMLO*(6LO>l>/9*3  002500 

LIHHI«(NHX/k)*l  002910 

IF  (NPTS.CO.flt  CO  TO  00  0C292C 

DO  90  JaNLO.JHXFO  CC293C 

VF  jVM*FLO»T(J-1»*OV  002990 

90  •' INT  900 » J,FF(J),VP  OC2990 

00  90  JaJJPO.NHI  002900 

XFa  VFT*FLOIT  ( J-l>  *OV  002970 

90  .'«INT  900,  J,FF< J> , VP  002900 

60  roNTINVC  uL  2990 

NLIH*NHI«NLO*l  002909 

VlPayFT*FLO»T(NLO-l»»OV  CC  25  10 

V2P*VFT*FLO»T(NHI-l>«OV  Ok  2929 

C VI*  IS  FINjrr  FKFO  OF  FONCLV  0C  2539 

C V ?P  IS  L»ST  F*CO  OF  PONE CV  0O259C 

* MNITC  (KFXLFl  VlP,V2P,0VP,NkI9  0C295C 

t HIIITE  <KFI«.E)  (FFCJl  ,>NL0,NHI)  012569 

PkFFEt  OUT  (KFILE.1I  (VtP.NLIN)  0C257C 

IF  «UNIT(KFILE1  .CO.  0 > STOP  OP250C 

9UFF£*  OUT  (KFILC.H  (FF  (Hl.01 , FF(  NHI)  ) 00259C 

IF  (UNITtKFTLE)  .CO.  0 > STOP  00 26 OC 

vft*  «,7*floot(nlipf-i>»ov  atibit 

IF  (1ST OP. CO. 1 1 CO  TO  190  0T2620 

JF«1  002630 

00  00  J*NIIHF,M»XF  0C2693 

FF(JF)«FF(JJ  C02650 

10  JFaJF.l  002660 

no  90  Jajr,H»XF  002670 

90  FF( JT«f . 002660 

NLC*MSHIFT*1  0C2690 

196  COLL  SECOND  (TIME)  002700 

TIHPNL*  TIHPNL  ♦TXHE“TINC0  002710 

WTU9H  OC  2720 

C 002730 

t 002790 

C 00279e 

900  F on  NOT  (110,29V  ,112.9, FI  2. 5)  0(2760 

CNO  002770 
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THIS  PAGE  IS  BEST  QUALITY  PRAGKCAfclil 

from  oopy  fuknishjsp  ro  ddq 


SUBROUTINE  MIRACL  76/76  OPT-1 


FTN  13/  16/TR 


1 


! 


10 


IS 


tl 


* s 


30 


3» 


60 


50 


33 


60 


63 


TO 


SURROUTINT  M IRACL  OCOIOO 

e oooiio 

C 000120 

000130 

• . 800160 

* • 0C01SC 

* CALCULATES  MONOCHROMATIC  ABSORPTION  COEFFICIENT  POP  SINGLE  LATER  * 000160 

* * 000110 

* * 001100 

• USES  HIRACC  ALGORITHM  FOR  THE  LCRENTZ  PROFILE  • 000130 

• * 000200 

00  0210 

C 000220 

C 000230 

COUPON  CNU(230>, 8(238), ALFAO(2SOI,£PP(23B7,MOL(230>,  000260 

C EF0PTNI236I »REC*LF|?30>  000250 

COMMON  FF  (36001  ,SF(900I  * VSF(223I  00C26C 

COMMON/PAIN/  KFILCtKPANEL.PO. TEMPO  000270 

COPPON/NEN/  NMOL <X ID (7) ,SECANT,PAVE ,TAVE« HHOLIO  (7)  *MK(  ?> , 000260 

C DV, VI, V2 ,MLATER  JC  0290 

COMMON  /XSUB/LlMIN,IL0,IMT,V607,VT0P,VFt,IE0F, IPANEL  ,ISTOP,IOATA  00C33C 
COMMON  /MAX/  NPF,OXF ,NF  ,NWS,OXS,NS,N WVS,DX VS , NVS »NFHAX»NSHAX»  NVSH  SC031C 
1 AX  0C072D 

COMMON  /SU91/  N AXF,MAXS  ,*)AXV$  ,NLIMF»NLIMS  , NL IPVS,  NLOtNPI,  OVS,  CWS  000330 
COMMON  /XPANEL/  VIP, V2P, OVP,NLIM,NSNIFT ,NPTS  000360 

COPPON/XTIME/TIPE,TlPROF,TIMCNV,TIHPNL  OP  0330 

DIMENSION  ALFCOP(T) , ALFOl ( 7) , U (?) , SCQR (7)  00CJ6D 

OIMENSION  XF (231) ,XS ( 2S1) ,XVS (231 > 000370 

C OC 1390 

OATA  NPF  / 6/,  PXF  / 0,02/,  NF  /201/,  NFMAX  /251/  uiuSOO 

OATA  NWS  /16/,  OXS  / 0.00/,  NS  /201/,  NSNAX  /25\/  OODoCO 

DATA  NWVS/66/,  OXVS/  0.32/,  NVS/201/,  NVSMAX/251/  OiOLK 

C IRU620 

DATA  IFNTER/0/,  LIPIN/230/,  NSMIFT/32/;  NUN."^  , 01/ , NPTS/  0/  000630 

c 000440 

DATA  SUBIO  /10M  MIPACL  / 0U45, 

C 300460 

IF  (IEPTEP.NE.OI  GO  TO  10  0(  0470 

IEMTE  R**l  OOOASj 

NLIMS»( NLIPC/6I  ♦ 1 uOUDC 

NLTMVS*  (NLINS/4)  ♦!  it  .5  30 

C NOTE  (DXVS/OXF)  IS  16  AND  (OXVS/OXS)  IS  6 000510 

NtlOUNO-FLPAT(NNVS>*(OXVS/OXF>  Of  0*20 

MAXF«NLIMF4N80UN0  C10530 

MAXS-IPAXF/6)  ♦!  000560 

MAXVS«(MA*S/6»6l  000330 

CALL  SH  APEK  XF, XS, X VS)  000360 

CALL  NOLECCi, HHOLIO, NMOL, TEMPO ,TAVE ,P0  * PAVE, SCOP, ALFCOR, ALFOl I 000570 

10  PRINT  300  0005SC 

PRINT  3TP,  SUMJp  01 0590 

tinpof»timcnv*ti:npni»o,  sco&j: 

PRINT  906,  (XTOI T>  ,1*1,7)  000610 

REMIND  XF1LE  010620 

IEOF*«  00 06 SC 

PRINT  *10,  KFILE  000660 

PRINT  919,  SECANT  000630 

PRINT  923,  *AVE,TA VF  CC0660 

PRINT  930,  CNMOL ID( Ml  ,NR(M)  ,N»1 , NMOL)  000670 

PRINT  960,  DV,V1,V2  iilObSG 

OVP«OV  OC  Ob  90 

OVS* tOXS/OXF>*OV  00  or  00 

09VS*  (DX9S/DXF)  *OV  000710 

BOUNO»FLOAT (N90UN0) *OV/2 • 060710 

BOUNVS-MUNO/2.  300739 

PRINT  962,  OOUWS  000760 

ALFNAX*  OOUNO /FLOAT  (NWVS)  000730 

NLO*NBOL.IO«1  000760 

NMI«MLIMP6NSMir*-l  000770 

00  30  1*1, MAX*  OC07BO 

30  FFtI)»0,  30  0 791 

00  60  I«1,MAXS  OUR30 

60  SritlaO.  Ok  06 10 

00  7*  1*1  »PA XVS  000  620 

70  VSF(I»«P,  0(0630 
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*HIS  page  is  best  quality  mcncABLl 
IROIl  COPY  FURBISHED  TO  MQ  — — ' 


suokoutine  nimcl  79/79  opt«i 


UK  t.MWI 


n/u/ri 


NFITE  IXFILEI  ( XIO(I>  ,l«l,  77 , 9EC9NT,PAVf  ,TAVf  ,(NHQLI0(N>  ,"al  ,7)  , 

t i*mi>,i«i,h*oo,ov,vi,*7 
•uffp«  outckfili.ii  (xiocii,nlatex> 
if  (unitckfil7i  .10.01  stop  -or* 

xft»vi-  »,**ouno 

V9OTa«.«OUN0 

¥TCIta«*OOUN0 

XKT0a|.H*l*TlNP6 

««r  •i.nimtn 

XKTPACa(l./XKT}  • (1./XKT0) 

CALL  NOLfCCI, *HOLIO,NNOl,TENPO, TAKE, P6, PANE,  SCO*, ALPCOO.ALFtlll 
00  *0  Hal.NNQL 
U(i|aMK(NI*SCO*(NI*f SCANT 
ICNTa* 

SUNALPpI. 

NCNNCaf 

CONTINUE 

CONTINUE 

call  srcoNP  (tiheoi 

17  (ICOF.NE.O)  CO  Tf  100 
CALL  PCFILE(GNU,S,ALFAO,EPPiNOL> 
cm  WCONP  ITINEI 
T INPOFaTIMPOEdlNE-TINEO 

IMIE07.NE.»>  GO  TO  1 AO 

M00I7T  LINF  OAT  A FO*  TEMPERATURE , PRESSURE , ANt  COLUMN  CENSITV 

no  no  i»ko,tht 
n*  wol  < 1 ) 

EFOPTMCII  aj(  I)  *IMN7 
IF  <rrreVN»Il.iE.0.7  00  TO  |T0 
ICHI»|CNT ♦! 

AL7l«IL7A0  <11*  A17C0*(M 
IFIALFl.r.l.nVI  r,o  TC  110 

PRINT  999,  GNU(  Tl  ,$(  II  ,AL*A0(II  ,I17L,0V,M 

AL7L-CV 

NCHMCaNCMNCTl 

17  (AL7L.LF.ALFNAX)  60  TO  120 

PRINT  990,  6NU (I) ,5 <I),ALrA0(I> , ALFL , A17MAX,M 

ACMaAlFMAX 

NCNNGallCHNCM 

CONT INUf 

5UMAL7*  SIIMAL  F * AL  FL 
R!CAL7(I)al./ALFl 

FFOPTN(I7»rFOPTN<I)»EXP»-«‘PP<I»»XXTFACI»l»eC»L7(I>» 

* (l.-EXP(-GNU(II  /XXTIIFIt,-EXP(-r.NU»H/XKT07) 

CONTINUE 

IF  (NCHNG.6T.10PI  60  TO  160 
CONTINUE 

CALL  CONVFNL (GNU, EFOPTH, •ECALF»H3L*FF*SPtV3F , 

C X7,XS,XVS) 

tF(TPANn  ,f0.0)  GO  TO  100 

call  PANEL  ( FF, SF, V5F,*F ILf I 

GALL  SEPONO  (TIHEI 
X*ANEL»KPBNE171 
IF  (I5T0P.NC.1)  60  TO  1N0 
FNP  FlU  KFILE 
PONT  IMUE 

PRINT  999, TINE, Tt»*«>Or, TINCNV, T INPNl 
17  (ICNT.LE.01  GO  TO  170 
AV  »LF»5u» AL* /FLOAT  ( TCNT1 
POINT  99(1,  AXIL  F,  ICNT,  NCHW5 


0 0 00  9 0 
000660 
CO  I860 
000170 
100660 
OCOOOQ 
100900 
000010 
OCOUO 
000610 
000990 
000900 
000900 
000970 
0S09OC 
069990 
001000 
001019 
OUCKO 

ocioia 

001090 

6*1360 

601663 

onoro 

ooiooo 

CCU90 

3(1150 

nonio 

OvtlU 
001130 
311193 
001190 
0C116G 
001170 
5t 111  0 
301193 
0C173C 
Ot 12 1C 
Cl  1720 
0*1230 
CC129Q 
001290 
0(1261 
ocii7c 
0012*0 
0C129C 
cujw: 
061310 
or i3io 
001333 
OC 1390 
6C139C 
0C1363 
001370 
0(138  0 
351390 
001900 
001910 
011920 
001930 
001990 
0«  1990 
001960 
3(1970 
001980 
001910 
0619  00 
901*10 
001910 
06193C 
C i.  199  6 
001998 
OC  1960 
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THIS  PAGE  IS  BSST  QUALITY  PRACTICABLE 
mot!  COPY  PUBKIwHlU)  TO  W)C 


SHOWTIME  niracl  TWT* 


MK  k.ltttl 


03/H/T0 


ISO 

to* 

HO 

100 

iro 


c 

tro 

c 

c 

c 

000 

00* 

*10 

010 

too 

*20 

*90 

0*0 

0%t 

0*0 

0*0 

0*0 

0*0 

ooo 


MTU** 


format  i mot 

FORH*TtS*,?Al*t 

format  <*a  out*vt  fue  no.  • • i*> 

F0RNATIM  IICMT  Ml.lt 
FORMAT  <1*  (iHiO.Tt I 

format! ••  ooriitxit  ••  oit.o/  *e  t|NM«t  ••  rii.tt 

FORMAT ( 13H0  COtUMN  0CN*XTY<M0ttCUt.tS/CN**2t  , 

1 / /ton,**,  • • • lNtll.lt  t 

form*t<*i  d-ficm-d  • *Ptt.o/*o  ymcm-ii  ■ *fii.*/*o  vt<cn>it 


i *«.*» 

FORMAT ItfNI  B0UNn«0<CN>l)  ■ ,F*.*t 

rwott  <•  - — — — - 

rgiHf  t*  *4*4**444444***4  *444*4*4**4** ll.A ,t 1*. 3,3F10»*»I*t 


FORMAT  «lN#,10K,*TIt«»,ll*,*A*AO%*R,»CCN*Ot.UTION*,10K,*RANEt 
1 *rj*.9t 

FORMAT  <*0  AO  ERA SO  WIDTH  ■ *010.0, 

r • M>.  urns  • *110.  • NO.  CNANStJ  • »I10  t 

FORMAT  1 1M0 1 OK  .At  It 
TNP 


•/ 


OtlOTO 
OC  10  00 
0(1000 
it  1*00 
011*10 
otino 

0(1*30 
0(1**0 
0(1*00 
001**0 
001*00 
OC  1**0 
001**0 
Q01O00 
OtlOtO 

ot  tore 

*01030 

0*10*0 

0010*0 

ooioec 

001000 

00100C 

ottroo 

CUAnO 


SOS ROUTINE  C ONWFNL  0*/0*  OFT*! 


FTN  *• »**! A C3/l*/0* 


1 


* 


10 


10 


re 


r* 


3* 


3* 


*C 


*0 


SURROUTINt  CONVFNtlGNt'.EFnRTM.RECAtF ,MOL,FF  ,$0,VSF .KF.KS.K VS  t 
CCNHON/NCN/  NMOl  »XI 0<*t  , SECANT  ,R AOI.T AOf  , OHOLIDI  0>  ,NKC  0)  , 

C DO,  VI, *2, Mt AYER 

CON»ON  /VSU*/|TN1N,HO,IHI,V*OT  .VTftR  ,VFT  ,IEDF  ,I"*Ntl  ,1  STO*,I  DAT  A 
COMMON  FNXKF  UNO  ,OKF,NF,NNS,D KS,NS,NNVS  , OX  VS  » NVS,  NFNAX»NSHAX  , MVSN 
1 AT 

COMMON  rtUOl/  HIIK  0,  OAKS  , NXXVS  ,Nt  INF  , Nil  MS,  Nl  IMVS,  Nt  O.NHI  , DVS,  D VV3 
COMMON/  KTIME/TIHt,TIMFOF,TTMCNV,TIHRNL 

njM'NSICN  r.NUUI.  EFOFTH  (It,  RECAtFtll,  00  C 1 > , SMlt,  tSF(l) 
DIM»N*10N  MOUIt 

dimension  »Fin,  xs(t),  mm 

C At l SECOND  (TIME PI 
»ATV»*DVVS/0 »vs 
FVS*OtD*T  (NNVSl  /DOTS 
CON»*OVVS/OV-1. 

CONS*CVVS/PVS-l. 

IF  CUO.DT.  INIt  GO  TD  *0 
DO  *0  I*HO,  INI 
DE*TNl*  TFDDTNCI  l 
IFIDEFTNI.tt.Ott  60  TO  *0 
mCRI**ECAlF(It**ATVX 
»IMV*<6NU<Il-VFTt/0VVS 

irTS.»ys/iste»l 

JM*X*<IINT**FVS1*2. 

IF  t JMAX.LC. NAXVSt  60  TO  10 

U»*T.M 

GO  TD  *0 

20  JNIN*  tOINT-RFVSI  ♦ 

JMINF*<IINT*CONFt 

JMINJ»IOINT»CONSt 

0»S»(FtO*TWNTN-0t*01NTt*TStORt 

IS*IVS*<FtOAT<  JNINSt-IINT»CONSl*ISlORl 

IF* 1*34 | Ft OAT < JN INF !• II  NT* CONE  I* ISt ORE 

no  30  JJ* JNIN, JmAX 

JS*  JNIMS4JJ 

JPa  JMINF*  JJ 

IVS*IVS*I*tORE 

IS*  OS*OStDFt 

IF.tF.ISt  ORE 

T2V«*tRS<7V*l*l.* 

I0S*AOF <031*1.? 
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